SEMI-CLASSICAL GREEN KERNEL ASYMPTOTICS FOR 
THE DIRAC OPERATOR 



OLIVER MATTE AND CLAUDIA WARMT 



Abstract. We consider a semi-classical Dirac operator in d € INI spatial di- 
mensions with a smooth potential whose partial derivatives of any order are 
bounded by suitable constants. We prove that the distribution kernel of the 
inverse operator evaluated at two distinct points fulfilling a certain hypoth- 
esis can be represented as the product of an exponentially decaying factor 
involving an associated Agmon distance and some amplitude admitting a 
complete asymptotic expansion in powers of the semi-classical parameter. 
Moreover, we find an explicit formula for the leading term in that expan- 



1. Introduction and main results 

The free Dirac operator in d e N spatial dimensions is the matrix- valued partial 
differential operator given by 

d 

(f.f) D h>0 := a ■ (-ihV) + a := a k (—ih d Xk ) + a , he (0,1]. 

k=l 

The Dirac matrices ao, . . . , a<i appearing here are hermitian (d* x d^-matrices 
satisfying the Clifford algebra relations 

(1.2) {a k , a e } = 2 5 M 1 , k, £ = 0, 1, . . . , d . 

According to the representation theory of Clifford algebras such matrices exist 
and the minimal choice of their dimension e 2N is d* = 2^ d+1 '^. The special 
choice of the Dirac matrices is immaterial for our purposes; only the relations 
(If ,2p are used explicitly below. It is well-known that, as an operator acting 
in the Hilbert space L 2 (R d , C d *), D hfi is essentially self-adjoint on C^(R d , C d ") 
and self-adjoint on H 1 (R d , C d *). Its Fourier transform can be easily diagonalized 
(compare (12.61) and (12.71) below) revealing that its spectrum is purely absolutely 
continuous and given as 

(f .3) a(L\ ) = (Tjp h)0 ) = (-oo, -f] U [f , oo) . 
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Next, we add a smooth potential, V, to the free Dirac operator, 

(1.4) D hy := D hfi + Vl dt . 

We shall always assume that V has the following properties. 
Hypothesis 1.1. V G C°°(R d , R) and, for every multi-index a G Nq 7 

(1.5) sup \d™V(x)\ < oo. 

xeR d 

Moreover, there is some 5 G (0, 1) such that 

(1.6) - 1 + 5 ^ V{x) < -5, x G R d . 

In view of (11. 3p the previous hypothesis clearly implies that D h V is self-adjoint 
on i7 1 (IR d , C d *) and continuously invertible. In fact, its symbol, 

D v {x,0 := a-£ + a + V(x), (i,()eR M 
is uniformly elliptic in the sense that 

| det (D v (x, 0) | = (1 + |£ | 2 - V\x)) d * 12 > (26 - 6 2 ) d */ 2 > , (^)e R M 

Therefore, the inverse -D^y is given by some matrix- valued /i-pseudo-differential 
operator whose distribution kernel, R d x R d 3 (x,y) D^y(x,y), is smooth 
away from the diagonal. Our goal is to study the semi-classical asymptotics of 
this kernel, for fixed x ^ y. 

To formulate our main result we first introduce an associated Agmon dis- 
tance, d,A, on R d . It is the Riemannian distance corresponding to a metric 
conformally equivalent to the Euclidean one on R d , namely 

(1.7) G(x) := {1-V 2 {x))l d , xeR d . 
The Agmon distance is thus given as 

(1.8) d A (x,y):= inf / (q\G(q) qf' 2 , x,yeR d , 

where the infimum is taken over all piecewise smooth paths q : [0, b] — >• R d , for 
some b > 0, such that q(0) = y and q(b) = x. We also introduce an associated 
Hamilton function, 



(1.9) H(x,p) := -v 7 !- \P\ 2 -V(x), x,pe R d , \p\ < 1, 

and recall the following fact (see, e.g., j2j pp. 197]): If a smooth curve ( 7 ) : 
j _^ ^2d on a non _trivial interval / is a solution of the Hamiltonian equations 

such that 

(1.11) H(y{t)Mt)) = <>> 
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then 7 is a geodesic for the Agmon metric G. Let exp y : T y R d — > R d denote the 
exponential map at y G R d associated to the Riemannian metric G. We recall 
that two points x,y G R d are called conjugate to each other iff the derivative 
exp' y (v) is singular, where v G R d is chosen such that exp y (t> ) = x. In this article 
we shall restrict our attention to arguments of the Green kernel fulfilling the 
following hypothesis. 

Hypothesis 1.2. x*,?/* G R d , x* 7^ y*, and, up to reparametrization, there 
is a unique minimizing geodesic from y+ to x*. Moreover, x* and y* are not 
conjugate to each other. 



We recall that Hypothesis 11.21 is always fulfilled, for fixed y*, provided that x* 
is sufficiently close to y+. To state our main results we also introduce the - in 
general non-orthogonal - projections A defined by 



+ , : 1 „ 1 _ ct ■ ( + a 



1.12) A ± (C):=-1 ± -5(C), 5(C) :=-/==, CGC d ,|SCl<l 



compare Subsection 12. 1[ Here and henceforth we abbreviate Q 2 :— Cf + ' • - + Cj> 
for every £ G C d , and yf~- denotes the branch of the square root slit on the 
negative real axis satisfying ^t^/~- > 0. The following theorem presents the 
main result of this article in the case d ^ 2. 

Theorem 1.3. Let d ^ 2 and assume that V fulfills Hypothesis li.il and X*,?/* 
fulfill Hypothesis li.M Let ( 7 ) : [0, t] — > R 2d be a smooth curve solving (ll.lOp 
and satisfying (11. lip such that 7(0) = y± and j(t) = x*. Then, as h > tends 
to zero, 

n _ w si (i-v 2 (x*)) d -^a-v 2 (y*))^ (i + o(h))e- d ^y^ h 

ft det [exp^exp- 1 ^*))] 7 (2tt d A (x*, </*) A) 2 

(1.13) -f/(r)(-y(^))A + (^(0)), 

where U(t), t G [0,r], is a unitary matrix such that U solves the matrix-valued 
initial value problem 

27ie ierm abbreviated by 0(h) in (I1.13P admits a complete asymptotic expansion 
in powers of h. 

Proof. This theorem follows from (12. ip . (12. 2p . Proposition 16. 1[ and Lemma [6.41 
below. □ 

1 /2 

Remark 1.4. (i) The factor g(x,y) := det [exp^exp" 1 ^))] is familiar from 
the asymptotic expansion of the heat kernel associated to G, where it also 
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appears in the denominator of the leading coefficient. In particular, is it known 
to be symmetric, g(x,y) = g(y,x). 

(ii) It follows from Remark Wl\ that M(x^y±) := U{r) (-V(y±))A + (im{0)) = 
(-1/(x*))A + (ra(r)) a U(r) (-V(y±))A + (iza(0)). Using the latter formula we 
verify in the same remark that M(x±,y*)* = M(y+,x+) so that (11.131) has the 
correct symmetry property of the kernel of a matrix-valued self-adjoint opera- 
tor. 

(iii) In Appendix |A] we explain, in the case d = 3, the connection between 
the term U(t) (— V(y±))A + (im(0)) and the BMT equation for the Thomas pre- 
cession of a classical spin along a particle trajectory. The BMT equation is 
discussed in connection with the semi-classical analysis of the time evolution 
generated by D h y in [TJ [S] . O 

Next, we state our main result in the case d = 1, where we do not need any 
restriction on the entries x ^ y of the Green kernel. 



Theorem 1.5. Let x, y G R, x ^ y, and assume that V fulfills Hypothesis li.il 
with d — 1. Let (^) : [0,t] — > R 2d be a smooth curve solving (jl.lOp and 
satisfying (II. lip such that 7(0) = y and 7(7") = x. Then, as h > tends to 



1 

h 



(1 + 0{h)) exp ( - 








{l-V\x)) 


V4(l - 


^( y ))l/4 





(1.14) • ( cos(tf(r)) 1 - z sin(^(r)) (-F(y)) A + (iw(0)) , 

w/iere 

' io 2V( 7 (t)) 

T/ie term abbreviated by 0{h) in (11.141) admits a complete asymptotic expansion 
in powers of h. 

Proof. This theorem follows from (12.11) . (12. 2p . and Proposition 16.21 □ 



Remark 1.6. The choice of the sign of V in Hypothesis 11.11 is not important 
for our results. We could equally well consider smooth potentials V : R d — > R 
satisfying (II. 5p and 0<5^U^1 — 5. In fact, this immediately follows from 
the following trivial observation: If ao, ■ ■ ■ , cid are Dirac matrices and Dh.v 
is defined as in ( II. ip and ( II. 4p with some positive V, then Dhy = —Dh-v, 
where D^-v '■= oc ■ (—ihV) + 5o — V. Here the new Dirac matrices ctj := —(Xj, 
j = 0, . . . , d, again satisfy (II. 2p and, hence, Theorems 11.31 and 1 1 . 5 1 are applicable 
to Dh-v- In doing so we first observe that, for two given points £*,?/*, the 
validity of Hypothesis 11.21 does not depend on the sign of V since the Agmon 
metric G = (1 — V 2 ) 1 depends only on V 2 . Moreover, the expression in the 
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first line of the right hand side of f ll . 13j) . which we denote by A(x+, y+), does not 
depend on the sign of V either, since the Agmon distance and the exponential 
map are defined by means of G. We have, however, to introduce a new Hamilton 
function, 

H(x,p) := -v 7 ! - M 2 + V{x) , xeR d ,\p\<l. 
Let (!) : [0,r] be a Hamiltonian trajectory solving (ll.lOp and (11.111) with H 
replaced by H such that 7(0) = y± and j(t) = x+. Then the last line of (11.131) 
has to be changed as follows. Since -D^y = —D^ l _ v we obtain 

(1.15) D^ v (x*,y*) = -A(ar*,y,) U(r) V{y*)h7 (iw($)) , 

where the projection A~ is again defined with the original atj and U(t), t E [0, t\, 
is a unitary matrix such that U solves the matrix-valued initial value problem 

By Remark 11.4( h) we may multiply the right hand side of ( ll.ISp from the left 
with V(x+) A^(iw(r)) (— q;q). Similar replacements have to be made in Formula 
(11.141) for the one-dimensional case. O 

Example 1.7. Assume that V — E is some constant spectral parameter E E 
(— 1, 1). Then we can compute the Green kernel of Dh,E — a • (—ihV) + «o + E 
by means of the Fourier transform and find the well-known expression 

°^y) = {l^wr 12 - v\lh) 

(l - E 2 )^ 1 /r\W 



(2ir)* 2 h d 



(116) + ( ao - £ + iftW 2-D^) ^ ( ^g r/A) }. 

where we abbreviate r := |zc — y| and r := x — y in the second line. For 
large p, the Bessel function of the second kind, K d / 2 , behaves asymptotically 
as K d/2 {p) = {>K/2pyl*e-o{l + 0(1/ p)) and K' d/2 (p) = -(tt^p) 1 ^ c -p (1 + 
0(1/ p)). Moreover, it is clear that \/l — E 2 is the constant momentum of 
the Hamiltonian trajectory running from 7/ to x in the level set {p 2 = 1 — E 2 } 
and we readily verify that d^(x,y) = yl — E 2 \x — y\, (1 — V 2 (x)) l ^(l — 
V 2 (y)) l / 4: = \/l — E 2 , exp' y = 1, and [7=1. Consequently, the leading asymp- 
totics in (I1.16P agrees with the value predicted by Theorems 11.31 and 11.51 O 

An asymptotic expansion analogous to (11.131) has been derived earlier in [5] for 
a certain class of /i-pseudo-differential operators whose symbols are periodic in 
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the momentum variables. In the general case encountered in [5] the Agmon 
metric is replaced by a suitable Finsler metric. This is due to the fact that the 
figuratrix at x G R d , 

(1.17) f x := {peR d : H(x,p) = 0}, 

which is well-defined due to (11. 6p and just a sphere with radius a/1 — V 2 (x), is 
replaced by the boundary of some more general strictly convex body in more 
general situations. We also remark that the exponential decay of eigenfunctions 
and the semi-classical tunneling effect for the Dirac operator is studied by 
means of the Agmon metric in [TP] , 

We briefly outline the strategy of our proofs and the organization of this 
article. The first step is to conjugate the Dirac operator with exponential 
weights e v ' h where cp is essentially given as the Agmon distance to compare 
(I2.ip - (I2.3I) below. The choice of ip is explained more precisely in Section [2j Its 
construction is the same as in [5] and we shall not repeat the details of the 
proofs in Section [2j The distribution kernel of a parametrix of the conjugated 
Dirac operator yields the prefactor in front of the exponential in (11.131) and 
(11.141) . The symbol of the conjugated Dirac operator, whose properties are 
also discussed in Section [21 is given by a non-hermitian matrix having two 
different (d*/2)-fold degenerate eigenvalues, one with a non-negative real part 
and another one with a strictly negative real part. Only the part corresponding 
to the eigenvalue with non-negative real part contributes to the asymptotics of 
the distribution kernel. To obtain the asymptotics we first construct, roughly 
speaking, a parametrix for a "heat equation" (backwards in time for the part 
of the symbol belonging to the eigenvalue in the left complex half-plane) by 
means of a WKB construction. Since the eigenvalues are complex we use a 
Fourier integral operator with complex-valued phase function as an ansatz for 
the parametrix and work with almost analytic extensions. In order to solve 
the associated complex time dependent Hamilton- Jacobi equation we adapt the 
constructions of [HE]. In Section [3] we provide a self-contained discussion of the 
time dependent Hamilton- Jacobi equation that proceeds along the lines of [7] 
and provides some alternative arguments to control the derivatives of certain 
error terms and implicit functions. To solve the transport equations in our 
WKB construction we employ a strategy based on the Clifford algebra structure 
we learned from [TXj . This strategy is adapted to our setting in Section HI 
(WKB constructions for the usual Dirac equation which also apply to non- 
scalar potentials can be found in [31 IB]-) in Section we construct a parametrix 
for the conjugated Dirac operator by integrating the parametrix for the "heat 
equation" with respect to the time variable and adding a term accounting for 
the part of its symbol left out in the WKB construction. Finally, in Section [6] we 
compute the asymptotics of e^^ 11 D^y(x, y) e~^ y ^ h by means of a stationary 
phase expansion in the time variable and the momentum variables of the Fourier 
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integral operator. The main text is followed by an appendix where the BMT 
equation for Thomas precession is related to our results. 

2. Construction of a weight function 

2.1. Eigenvalues and eigenprojections of the symbol of the conjugated 
Dirac operator. Let us describe the first step in the derivation of the asymp- 
totics (11.13p and (11. 14p . We fix two distinct points, x* and y*, in R d fulfilling 
Hypothesis 11.21 and seek for some bounded weight function ip 6 C°°(R d , R) 
satisfying 

(2.1) cp(x*) - (p(y*) = cJa(z*,2/*) • 

Since (p is bounded and smooth it is then clear that 

(2-2) D;' v (x*, to) = e~^ h D^Jx*, y*) e^'\ 

where R d x R d 3 (x,y) h- >■ D h y ip (x,y) denotes the distribution kernel of the 
inverse of the conjugated Dirac operator 

D h>v , v := e^ h D hy e-^ h 

(2.3) = ot ■ (-ihV + iVp) + a + V l dt . 

The Green kernel asymptotics of Dh,v,tp thus yield the prefactor in front of the 
exponential in (11. 13)) and (11. 14j) provided that p is chosen in the right way. To 
motivate the partial differential equation determining p we observe that the - 
in general non-hermitian - matrix 

(2.4) D v (x,() := a-( + a + V(x)l, (x,()6R rf xC i , 
has two (d*/2)-fold degenerate complex eigenvalues, namely 

(2.5) X ± (x, C) := ± VT+C + V(x) , (x, C) € R d x C d , \%(\ < 1 . 

The eigenprojections corresponding to the eigenvalues in (12. 5p are given by 
(11.121) . In fact, a straightforward exercise using ( II. 2p . which implies (ex. ■ () 2 = 
C 2 1 and, hence, (ex ■ ( + a ) 2 = (( 2 + 1)1 reveals that, for ( e C d , \%(\ < 1, 

(2.6) A + (C) + A-(C) = 1 , S(C) 2 = 1 , A^C) 2 = A ± (C) , 

(2.7) D v (x,C)A ± (C) = AiCOA^C). 

Using cos 2 (#/2) = (l + cos(0))/2 ^ cos(9), 9 e (-tt/2, tt/2), we further observe 
for later reference that 



cos(#/2) ^ y/\z\ cos(#) = Vfti, z = \z\ e ie G C , %tz > . 

In particular, 



(2.8) & VI + C 2 > VI + (^C) 2 - (^C) 2 , CeC d ,|3C|<i 
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2.2. Agmon's distance as an optimal weight function. In this subsection 
we treat the partial differential equation determining the weight function ip. 
This equation is the eikonal equation corresponding to the Hamilton function 
introduced in ( II. 9p which is related to the eigenvalue A + as 

(2.9) H(x,p) = -v 7 ! -P 2 - V(x) = -\+{x,ip), x G R d , \p\ < 1 . 
We start with an elementary proposition covering the one- dimensional case. 

Proposition 2.1. Let d = 1 and assume that V fulfills Hypothesis Then 
the following assertions hold true: 

(i) For y G R, the unique solution in C 1 ([R, R) of the initial value problem 
H{x,(j)'(x)) = -y/l - <p'(x) 2 - V(x) = 0, xeR, ±0'>O, <f>{y) = 0, 

is given by the smooth function 

4>{x) = ± I y/l - V 2 (t) dt , xeR. 
Jy 

We have <f>(x) = ±d,A.(x, y), for x ^ y, and <f>(x) = +d^(x, y), for x < y. 

(ii) Given x+,y* G R, x+ ^ y±, we find some compact interval, K C R, such 
that x±, y+ G K and some p G C°°(R, R) such that <p(x) — (p(y) = d^(x, y), for 
all x,y G K with sgn(x — y) — sgn(:r* —y*), p is constant near ±oo ; and 

(2.10) H{x,p'{x)) < 0, xeR, and H (x, <p'(x)) = x G K . 

Proof, (i): Since -1+5 < V < -5, every solution <p G C\R, R) oiH(x, <j/(x)) = 
0, x G R, satisfies either <fi' > or 0' < on R. Thus, H (x, 0') = is 
equivalent to either <f> = yl — V 2 or <f> = — yl — V 2 and the first assertion 
is evident. Since the expression J( q \ G(q) q ) 1//2 is invariant with respect to 
reparametrizations of the path q, we may plug in q(t) — y + 1, t G [0,x — y], 
for y ^ x, or q{t) — y — t, t G [0,y — x], for y > x, to verify that <f>(x) = 
±sgn(x-y) d A (x,y). 

(ii): We choose some compact interval K with K 3 x+,y+ and pick some 
9 G C °°(R, [0, 1]) su ch that 9( t) = 1 if and only if t G K . Then we define 
(p(x) := ± f^8(t) a/1 — V 2 {t) dt, x G K , where we choose the +- sign if and 
only if x± > y±. By Part (i) p> satisfies H(x, p>') = on K and, for x K , we 
have p\x) 2 = 6(x) 2 {l - V 2 (x)) < 1 - V 2 {x), that is, H{x, p'(x)) < 0. □ 

To discuss the multi-dimensional case we denote the flow of the Hamiltonian 
vector field corresponding to H by $ = (X, P) : £>($) ->■ R 2d , so that X>($) = 
{(t,x,p) G R x R d x B x : t G I m3 ^{x,p)} and 
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Here I max (x,p) denotes the maximal interval of existence for the initial value 
problem p = (V p H(p), -V x H(p)), p(0) = (x,p), and B x = {p G R d : p 2 < 1}. 

In order to recall a result from |5j we remark that, if Hypothesis 11.21 is 
satisfied, there is, up to reparametrization, again only one minimizing geodesic 
running in the opposite direction from to y*. Moreover, we can prolong the 
geodesic from y* to x* or from to y± a little bit such that it still remains 
minimizing. Now, we are prepared to recall the following special case of [3J 
Proposition 4.5]: 

Proposition 2.2. Let d ^ 2 and assume that V fulfills Hypothesis \1.1\ and x* 
and y* fulfill Hypothesis \1.2\ Then there exist a point, y$, on the prolongation 
of the geodesic from x+ to y+, a compact neighborhood, K , of the geodesic 
segment from y+ to x+, some open set, W C T yo R d = R d , which is star-shaped 
with respect to zero, and some bounded function, (p G C°°(R d , R), with bounded 
partial derivatives of any order such that the following holds true: 

(i) For all x G R d , we have \Wtp{x)\ < 1 and 

(2.12) H(x, Vip(x)) < 0, and H(x, Vip(x)) = <=> x G K . 

(ii) ip(x) — <p(y±) = d^x^y*), for all x on the geodesic segment from y+ to 
x* . 

(Hi) exp yo \ w e C\W, R d ) n C°°{W \ {0}, R d ) is infective on If and 

K C ex Vyo (W) \ {y } . 

(iv) For every x G K , there is a unique pair (r,po) G (0, oo) x fi yo such 
that the projection of [0,r] 3 t i— > $(t,yo,Po) onto R d is a minimizing 
geodesic from y to x. We have 

(2.13) Hr,y ,p ) = (x,V<p(x)) = (X(r, y ,p ), V<p(X(r, y ,Po))) . 

Proof. We proved this proposition in (5] assuming that H G C°°(R 2d , R) and, 
for all x G R d , the function H (x, ■ ) : R d — > R is strictly convex, even, and 
H(x, 0) < 0. On account of (11.61) we can easily modify H such that these 
conditions are satisfied. To this end we first restrict H given by (12.91) to the set 
R d x {p G R d : \p\ ^ 1 — S 2 /2} and then we pick an arbitrary smooth extension 
to R 2d of this restriction fulfilling the condition imposed in [5]. The assertions 
of the present proposition do not depend on the choice of the latter extension. 
Another condition required in [5] is that 

(2.14) inf { F(x, v) : x G R d , v G S^ 1 } > , 

where F : R 2d — > R is the Finsler structure given by F(x, v) := (v\ G(x) v ) 1//2 = 
a/1 — V 2 (x) \v |, x,v G R d . Of course, (12.141) is a trivial consequence of ( II. 6p . 

□ 
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2.3. Consequences for the symbol of the conjugated Dirac operator. 

In this subsection we collect some important properties of the eigenvalues of 
the symbol 

(2.15) D Vjtp (x, := D v (x, f + iV<p(x)) , (x, £) G R 2d , 



always assuming that V fulfills Hypothesis II. 1| x* and y* fulfill Hypothesis ll.2| 
and that <p is the function provided by Proposition 12. 1| for d — 1, or Proposi- 
tion I2.2| for d ^ 2. We introduce the complex-valued symbols 

(2.16) o±(z,0 := TiA ± (x,e + zV^(x)) 

= -i^TT (e + ^(x)) 2 =f , (x,0 G R 2d . 

Notice that, since \ V(p\ < 1, a + and a_ are well-defined, complex- valued smooth 
functions on R 2d . In the next two lemmata we collect some basic properties of 
a± which are used in the sequel. Henceforth, we abbreviate (a±)"^ '■= d^V x a±, 
Hp X := dxVpH, etc. 

Lemma 2.3. For all x, £ G R d , 

(2.17) Sa+(a,0 < 0, 

(2.18) 3a+(x,0 = ^ (x, £) E K x {0} . 
Moreover, we have, for all x G R d , 

(2.19) a+(x,0) = i#(x,V^(x)), 

(2.20) V e a+(x,0) = V p H(x,Vip(x)), 

(2.21) (a + )l(x, 0) = fl^s, V<p(x)) + H'^{x, V<p{z)) <p"(x) , 

(2.22) (a+)jk(z, 0) = -iiZ£(z, Vp(z)) . 
In particular, for all x G i^o ; 

(2.23) a+(x,0) = 0, V x a+(x,0) = 0, (a+)^(x, 0) = , 

(2.24) V^a + (x,0) = -Vcp(x)/V(x) ^0. 

Proof. By virtue of (12. 8p and Proposition I2.2H ) we obtain, for x, £ G R d , 



3a + (x,0 = -3?v / l-(V^(x)) 2 + e + 2i(e|V<^(x))-^(x) 



«C -V / l-(V<^(a;)) 2 + £ 2 -\/(a;) ^ 0. 

Moreover, by (12.121) and the previous inequalities, 5$a+(x,£) = if and only if 
£ = and x G Kq, which yields (12.171) and (12.181) . Furthermore, V p H(x,p) = 
—V p X + (x,ip) = —iV^X+(x,ip), for x G R d , \p\ < 1, whence V p H(x, Vip(x)) = 
V c a + (x,0), for every x G R d , which is <^M) . ( 12T2Tj) follows from (12T2Uj) . Next, 
Hp p (x,p) = —id p V^\ + (x,ip) = (X + )'^(x,ip), x G R d , |p| < 1, which implies 
(I2.22p . All remaining identities are obvious. □ 
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Lemma 2.4. For all x, £ G R d , 

(2.25) ^ -25 < 0. 



Proof. Since if (a;, Vy?) ^ on R d the inequality (I2.8P implies 



3fa_(a;,0 = -^V 1 = (V<^(a;)) 2 + £ 2 + 2z ( £ 1 V^(x) ) + V(x) 



< -V 7 !- (Vv?(x)) 2 + V(x) ^ 2V(z) ^ -2 5, 



for all x,£ e R d . 



□ 



3. The time-dependent complex Hamilton-Jacobi equation 
In this section we solve the complex Hamilton-Jacobi equation for ip±(t,x,rj), 



To this end we proceed along the lines of the constructions in [7|. Instead 
of solving the problem first for symbols which are homogeneous of degree one 
as in [7] and then using a standard reduction to that case (see, e.g., [9]) we 
carry through all constructions for general symbols as in [3]; see also [S] for 
symbols which are periodic in the momentum variable. Since a± is complex- 
valued we can only hope to solve (13.11) up error terms 0((Qtp±) N ). These 
do, however, not any harm in the WKB construction as we shall see in the 
proof of Proposition 14. 8} see, in particular, (I4.4ip . In order to control the 
derivatives of these error terms and of certain implicit functions appearing in 
the constructions we give some arguments alternative to those in [7]. The final 
result of this section is Corollary 13. 101 where (I3.ip fe fl3.2p is solved. 

In the whole Section [3] we always assume that V fulfills Hypothesis 11.11 and 
x*, y* fulfill Hypothesis 11.21 ip is the function provided by Proposition 12.11 for 
d = 1, or Proposition 12.21 for d ^ 2. 

3.1. Estimates on the Hamiltonian and contact flows. From now on the 
symbols x, y, £, and rj will denote complex variables in C d . We set p = (x,£) 
and d Xj = (d^ Xj - id% x .)/2, <).,._ = (d^ Xj + id^ Xj )/2, V, = (V : ,,v - iV% x )/2, 
Vf = (Vsfcj; + iVc$ x )/2, . . . , and analogously for complex variables other than 
x. In the rest of this article we further extend (p and V almost analytically to 
smooth functions defined on C d - again denoted by the symbols (p and V - so 
that, for every compact subset K C C d and all N e N, a 6 N 2 ^, we find some 
CN,K,a £ (0, oo) such that 



(3.1) 
(3.2) 



^± + a ± (x,V^±) = 0({Sh/>±) N ) , iVGN, 
1/)±(0,X,T)) = (r)\x), Qtp± ^ . 



\ d (9tx,Zx)Vx<f(x)\ ^ C N>K ,a\3x 
^,3x)VxV r (x)| < C NtKta \$fX 



x e K , 



x e K . 



ii 



Then we find some open neighborhood, ft C C , of R such that the symbols 
a±(x, = -iy/l + (I + iV V {x)Y t iV(x) , (i,()gO, 

are well-defined and almost analytic on ft. For every compact subset K C ft 
and all iV G N, a G Ng d , we find some Cn,k,ci G (0, oo) such that 

(3-3) |9(V,%) V P a ±(P)l < C N,K,a \$*pf, P^K. 

(In fact, V^a± = on ft.) In this subsection we proceed along the lines of [U [7] 
to obtain estimates on the Hamiltonian and contact flows associated with a±. 
On ft we introduce the Hamiltonian vector fields 

jT a± := (V s a±|V x )-(V x a ± |V e ), 

and the elementary actions 

*4(P) := ( V 5 a±(p) | O - o±(p) , p = (x, G ft . 

Here and henceforth ( • | •• ) denotes the extension of the Euclidean scalar prod- 
uct to a bilinear form on C d . For later reference we infer from (I2.23P that 

(3.4) *f+(x,0) = 0, xeK . 

We further add an extra variable, s G C, to (x, £) G ft which parameterizes the 
action and define the contact fields 

J^ ± := -tf± d s + ^ a± = (a± - ( V s a± | £ )) 0, + JC± on C x ft . 

Finally, we introduce the real partial differential operators 

Notice that, since c d z + cd? = (JRc) d^ z + (3c) d<$ z , the vector in C 2d corre- 
sponding to J^ a± under the identification d^ Xj a,-, d^ x . <-> iej, dg%. <-> e^+j, 
j = 1, . . . , d, where (ei, . . . , e2d) is the canonical basis of C 2d , is 

just 

(3,) ^ „ (_^) . 

We denote the flow of as /tf = (Q ± ,S ± ) : D(k ± ) -»■ C M , so that 

£>(k±) = {(t 7 y,r)) G R x ft : (G </m ax (y, ??)}, where J^Av^) denotes the 
maximal interval of existence for the 4<i-dimensional real initial value problem 
P = J%l±(p), P(0) = d/,77), and 

- (-^£3>) 
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on P(k ± ). The flow of Jt a± is then given by (? ± ,k ± ) : C x T'(k ± ) C 1+2d , 
where 

<^{s,t,y,ri) := <£(s,y,rj) := s- / £/ ± (nf(y,ri))dr, s G C , (t, y, 77) G £>(«;*). 

Lemma 3.1. Lei xo € -Ko ond / C /max(£o, V<^(xo)) some interval such that 
X(t, xq, Vip(xo)) G Ko; / or t G /. T/ien / C </m ax (£0) 0) and 

Q + (t,x ,0) = X(t,x , V<p(x )) , E + (t,x ,0) = 0, tel. 

(Recall the notation introduced above (12.111) . ) 

Proof. By <^M> and (T2T231 we have JC + = V p F(x, V<p) ■ Vjh, on AT x {0}. 
Moreover, <9iX 4 (:r , V<p(2 )) = V p H {$ t (x , Vy?(s ))), and $ 4 (x , V<p(:ro)) = 
(X t , V^(X t ))(ar , V^(x )), t G J, by □ 

Lemma 3.2. Let r > and assume that p : [0, r] — > R 2d is a real integral curve 
ofJ^ + with p(0) G K x {0} = {3a+ = 0}. T/ien p([0,r]) C K x {0}. 

Proof. By assumption S(V^a + , — V x a + )(p(t)) = ^ £yp(£) = 0. Since a + fulfills 
the Cauchy-Riemann differential equations on the real domain, it follows that 
(^0+)^ (p(t)) = 5s(a + )' (p(t)) = 0. Hence, the derivative of (3a + )f| R 2d vanishes 
along p, thus 9ta+(/?(0)) = implies p([0,r]) C {Sa + = 0}. Using (I2.18p we 
conclude p([0,r]) C AT x {0}. □ 

In what follows we consider the trajectories of J^ a± emanating from the planes 

£ (r?) := {(-My,v),y,v y Mv,v)) ■ yec d }, 

where ipo(y,v) — (v\y)i V £ C d , 77 G R d , so that V y il>o(y, rj) = 77. The reason 
why we restrict our attention to real 77 is that in this case ipo trivially fulfills 
the inequality 

(3.7) 9^ (y^)-(%|»V^o(y^)) > -0(|^(y,r/)| 3 ), yGC d , 

which is used to derive the estimates of Lemma 13.31 below. (We could equally 
well consider more general ipo satisfying (13. 7p locally on compact subsets in 
Lemma EH) We define 6 : C 1+2d -> R by 

(3.8) 6(s,x,0 := -Qfs- = -3f(s + (s,x,OeC 1+2d , 

which corresponds to the function — ( | ) considered in [7] where the sym- 
bol is homogeneous of degree one in £. The following lemma is a slight modifi- 
cation of [T) Proposition 3.1]. For the convenience of the reader we present its 
proof in Appendix [BJ 
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Lemma 3.3. Let 2/0,770 G K d - In the minus-case we set r = 0. In the plus-case 
we pick some r G J^JdJo-, t]o), t ^ 0. Ifr > we assume that Kt~(yo, fjo) is real, 
for all t G [0, r], and ^sa + (y , r/ ) =0. Then there exist e > 0, C G (0, 00), and 
some neighborhood & C C 1+2d of (— ipo(yo, T]o), yo, Vo) such that the following 
inequalities are satisfied on £0(770) H /or a//0^r^t^r + £ and /i G [0, 1] ; 

(3.9) ) ^ - / -9fo±(^)du-C|9fK^| 3 , 

2 Jo 

(3.10) |^ ± | 2 + 6(q t J ^ ± )>^{|9f«±| 2 + 6(^,K±) + y' -3a±(3&c±)du}, 

(3.11) |^±| 2 < cd^ + e^,^)). 

The constants e and C can be chosen uniform when rjo varies in some compact 
set. 

In particular, if (s,p) G £o(?7o) H & and nf{p) and qf{s,p) are both real, for 
some t G [0, r + e\, then nf{p) is real for all r G [0, t] . 

We recall that the Hamilton matrix, F a± , of a± is given as (recall (13. 5p ) 

<-) r --(4 4)(S + (i^)ft> 

where 9 = {9 X , 6^) G C 2d and where the matrix on the right vanishes on the real 
domain R 2d C fi. Here and henceforth we write a'^ = |(c$ft£ + idQ^XVsfa; — 

iV<$ x ) a, etc. Let : C 2d -> C 2d denote multiplication with z. Then ( l3~B~j) 
implies that, for every compact K G Q and all iV G N, a G Ng d , there is some 
Cn,k,o £ (0, 00) such that 

(3-13) l\\d^ p) [y,¥ a ] 
on K, where a again is a + or a_. 

Corollary 3.4. Lei (yo,Vo), T , and e be as in Lemma fX3j set T := r + e, and 
/e/j X C C M 6e some sufficiently small compact neighborhood of (yo,rjo)- Then, 
for allN G N and a G N^ 4 " 1 , /3 G N^ 1 , t/iere exzsi C N)K ,T,a, C' NKTJ3 G (0, 00) 
swc/i i/iai 

(3-14) ||fl^ p) d^(p)|| < C w> « sup \%Kf(p)\ N , peK,te [0,71, 

se[o,t] 

and, /or (y,n) G LT n (C d x R d ) ; and t G [0,T] ; 
(3-15) 



a 



II 



as 

a" 




C N ,K,a \ ^SP 



N 
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Proof. Again we drop all ±-indices in this proof. If denotes multiplication 

by i, we have 

\\ d (t,sftp&p) d P K t(p)\\ = \\\ J i ^mp^'Ap) ]|| , p £ K. 

Moreover, we know that K' t (p), p G K, satisfies 4 «£(p) = F a (Kf(p)) «4(p), 
i G [0, T], «^(«J(p)) = k4 +s (p), £, t + s G [0, T], and n' (p) = 1. In particular, 
[j^ , ^o(p)] = 0. Since we have 

^X(p)] = [y, F a (K t (p))K' t ( P )] 

= F a (^(p))[^,^(p)] + [^,F a K(p))K(p), te [0,T], 
it thus follows from Duhamel's formula that 

[^X(p)] = / ^_ s (p)[^,F a (^(p))]<(p)rf S , tG[0,T], P 6K 

Using sup tg[0T] sup pgA - |K(p)|| < oo, and \J , F (/« s (p))] = OflQTc,^)!*), iV G 
N, we deduce that the following estimate is satisfied in the case a = 0, 

(3.16) \\[y,d^ p) <(p)]\\^C' ma sup \%K s (p)\ N , P EK, 

se[o,T] 

for some C' KNTa G (0, oo). If (13.161) holds true, for all multi-indices of length 
^ n G N , and |a| = n + 1, we write 

(3.17) j t \J, d^ p) K' t ] = ¥ a {K t ) \J, flg^ 9p) «a + t^> Fa(«t)] 5(V,3p)< 



Here we know that sup t6 [ 0T ] sup peK ||9^ ^ < oo, for ^ f3 ^ a, and 

II l^.^p,3p) F «(^)] II = 0(\%K t (p)\ N ) by ([313]). By the induction hypoth- 
esis || [J, 9jgJ> p )«t] II = <9(sup s6[0jT ] ^^(p)!^), for < /3 < a. Applying 
Duhamel's formula once more, using 9^ q \k' q ] = 0, we obtain (I3.16p . for 

a G Nq^, |a| = n + 1. Taking successively time derivatives of (I3.17P we can 
apply a bootstrap argument to include higher order time derivatives and to get 

(3.18) \\[J,df mp ,z p) K'M\\ SU P IWp)^ P^K, 

se[o,T] 

for some C'^ N ^ a G (0, oo) and all a G N^ d+1 . Combining fl3TT8|) with (EOT]) 
we arrive at the asserted estimates (13.141) and (I3.15p . □ 
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3.2. Approximate solution of the complex Hamilton-Jacobi equation. 

We recall the notation (I3.6P and set 

^+ := {(UO)GRjxif x{0} : Q+(y, 0) € K , f G [0, t] } , 
® + := { (t,x,0) G R+ x A x {0} : Qt f (x,0) EK ,t'e [0,t] } , 
g+ ■= ({0} x R 2d ) U <f + := ({0} x R 2d ) U 

^ + := {(t,g t + (y,0),0,y,0) : (t, y, 0) G § + } U {(0, x, rj, x, rj) : x, V eR d }. 
We can represent as the graph of the function (k + ,g + ) : S' + — » R 2d , where 

(3.19) fc + (*,z,0) = g + (-t,x,0) = X(-t,x, V^(x)), (t,x,0)G^ + , 

(3.20) k + {0,x,r]) = x, (x,rj) G R 2d , 
and 

(3.21) (? + (t,x,0) = 0, (t,x,0)G^ + , 9 + (0,x,rj) = rj , (x, rj) G R 2d . 
In the minus-case we simply set 

g~ := S~ := {0} x R 2d , &~ := {(0, x, r], x, rj) : x,rj G R d } . 
Then is the graph of (k~,g~) : c?~ — )■ R 2d , where 

(3.22) fc-(0,x,77) = x, g-(0,x,r)) = r], (x,r?)GR 2d . 

The next lemma shows that we can extend fc ± and g 1 * 1 to smooth functions 
defined in a neighborhood of <^ ± which represent the canonical relations given 
by the flow of J4f a± as a graph in the vicinity of 

Lemma 3.5. There exist open neighborhoods, & ± of <S^ in Rq x C 2d and 
je± ofF* in R+ x C 2d x n, and k ± 1 g ± G C°° , C d ) , such that, for all 
{t,x,£,y,r ] )eJ? ± , 

(x,0 = (Q?(y,7)),Ef(y,rj)) & (y = k^fax,^) A £ = 9 ± (t,x,v) ) ■ 

Proof. We define 

(3.23) F ± (t,x,£,y,r}) := (x - Q ± (t, y, rj) , £ - E ± (t, y, rj)) , 

for all (x, £, t, y, rj) G C^xD^). In the following we regard F as a R 4d - valued 
function of 8d+l real variables. At t = we have F ± (0, x, £, y, 77) = (x—y, £—77) 
and it is trivial that (-^ 1± )(g^ ! cj|,SR2/,%)l*=o : — ^ * s invertible. This already 
proves the assertion in the minus-case. In the plus-case we have, for general 
(x,^t,y,r,)eC 2d xV(K + ), 

(^ + WeW*>*>^) = (i d &y,v) ■ 

\-U-2d I- /(%,%)/ 
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We know, however, that Q + (t, y,0) = X(t, y, V<p(y)), for (t,y,0) G where 
we use the notation introduced in the paragraph preceding (12.111) . In view of 
dSH we further have (3Q+)^(t, y, 0) = (KQ+)^(t, y, 0) = dy[X(t, y, V<p(y))] 

and (3fJQ + )' % (t,?/,0) = -(3Q+) , Rj ,(t, y, 0) = 0, for (t, y, 0) G 0+ Likewise we 

have E+(t,y,0) = and, hence, (E+)' { ^ y) (t, y, 0) = 0, for (t,j/,0) G 0+. It 
follows that 

(^WW*.<?*<'.».°>.<>.».0) = (£ -Mnt,y,Vj(y))\ 9 i,j 

for all (t,y,0) G Moreover, the matrix d y [X(£, y, V</?(y))] is invertible as 
a time zero fundamental matrix of some matrix- valued ODE. In fact, it holds 
X(0,y,V<p(y)) = y, thus d y [X(0,y,V(p(y))) = 1, and 

d t d y [X{t,y,Vcp{y))] = E{t,y)d y [X{t,y,Vcp{y))), 



:X(t,y,V<p(y)) ■ 



for every (t, y, 0) G where 

(3.24) B(t,y) := V<p(x)) + H'^(x, V<p(x)) </(x)) \ x __ 

Since can globally be represented as a graph we conclude that functions 
k + and g + with the properties stated in the assertion exist and are unique if 
the neighborhood £f + is chosen sufficiently small. □ 

In the remaining part of this subsection we show that a solution of (13. ip is 
given by the following formula well-known from classical mechanics, 

ip±(t,x,r}) := (Z± o K±)(t,x,rj) 

(3.25) =(k ± {t,x,7 ] )\r ] )+ [ £/ ± (Kf(k ± (t,x, V ),r ] ))dr, 

Jo 

for (t,x,rj) G where 

Z±{t,y,r}) := -^{-{y\r}),y,r]) = {y\rj)+ / £/±(k± (y , r])) dr , 



for (t,y,rj) G V^k^), and 

(3.26) K±{t,x,ri) := (t,k±(t,x,7)), V ) , (t,x, V )e& ± . 

To this end we further introduce the following canonical weights which are used 
to control the error terms, 

f± := |^K ± | 2 + e(-Z±,K ± ), oiiDf^), 

r± := T±oK± = |S(x,^ ± )| 2 -(^x|% ± )+^V±, on 

Here we used that K,f(k zi: (t,x,r]),r}) = (x, g ± {t 1 x, 77)), for all (t,x,rj) G 5f ± . 
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Lemma 3.6. There is an open neighborhood, jV± C IRq x C d x R d , of the closure 
of S ± such that, for every compact subset K C JV±, we find some Ck £ (0, oo) 
such that, for all (t, y,n) £ K and r £ [0, t], 

%Z ± (t,y, V ) > (^I^H^r?)-^ Za ± (MK s (y, V ))ds 

(3.27) -^|S(g ± ,S ± )(t,y,r / )| 3 , 

(3.28) ~ |S(g ± ,H ± )(r,y,r / )| 2 < r ± (r,y,tj) ^ C K f ± (t,y, V ) , 

Proof. We apply Lemma I3T3"| for every fixed r] £ R d , recalling that the constant 
C appearing there can be chosen uniform when rjo varies in a compact set. In 
the minus case we always choose r = in Lemma 13.31 In the plus case we 
choose r = 0, if (yo,Vo) ^ K x {0}. If, however, rjo = and yo £ K , 
then we choose r = max{t ^ : X(r, yo, Vip(yo)) £ K , r £ [0, £]}. Then all 
assumptions of Lemma l3.3l are satisfied because $$a+(yo, 0) = and ^K t (y , 0) = 
5s(X(t, yo, V^(no)), 0) = 0, t £ [0,r], by (12.181) and Lemma l3TT| respectively. 

□ 

First, we derive some estimates on the derivatives of the implicit functions k^ 
and g ± . To this end we put 

(3.29) ^±:=Kl\ji), 

where K± is given by (13.261) and jV± by Lemma 13.61 so that JV± C ^ ± is a 
neig hborhood of S in R+ x C d x R d . 

Lemma 3.7. Let k ± ,g ± £ C°°( < 3 ± ,C d ) fre i/ie implicit functions provided by 
Lemma 1 5*. 51 Then, for all compact subsets K C ^/f^ and all N E N and 
a £ NQ d+1 , there is some CN,K,a £ (0, oo) stzc/i ina£ ; /or all (t, x, rj) £ K , 

\\^t,^,^,^) ( y, k± ^){t i x } Ti)\\ ^ C NiK!a r±(t,x,r]) N . 

Proof. Dropping all ±-indices and using the notation ( I3.23|) we have 

where all derivatives of Q and S are evaluated at (t, k(t, x, T]),r)), for (t, x, rf) £ 
. We denote the above matrices as A, B, and C, so that A = B^ 1 C . Let ^ n 
represent multiplication with % on C d = R 2d , that is, J^ n = l n ® (° ~ ). Writing 
[J^, A] := JidA — AJ d , etc., we then have 

(3.30) \J, A] = [J, B- 1 ] C + B- 1 \J, C\ = B- 1 [B, J\ B' 1 C . 
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Taking derivatives of ( I3.30p we obtain, for a G No +3d , 

for some combinatorial constants c(/3,7, 8) G (0, oo). Here the commutator 

i d lmx^x, V ) B ' J \ has the form (p = (y> v)) 

5Z^KpSp)-^' ' (Polynomial in the partial derivatives of k) . 
We know from Corollary 13.41 and (I3.28P that 

II 

□ 

In order to show that the formula ( I3.25P defines a solution of (13.1 p we adapt 
a standard proof from classical mechanics and compare the differential of Z± 
with the pull-back under the map 

e±(t, y, rj) := (t, K±(y, V )) = (t, Q ± (t, y, 77), E ± (t, y, V )) , (t, y, 77) G V(^) , 
of the Cartan form, 

w± := £dx — a±(x,£) dt . 

u± is considered as a form on R x C 2d , so that dxj = d^txj + id^Xj, and we 
abbreviate £ dx := £1 dx\ H — • + ^ da^, etc. 

Lemma 3.8. On every compact subset K C ©(ft 1 * 1 ) swe/i i/iai |t| ^t^ on K 
we have, for all N G N and a G Ng d+1 ; 

( 3 - 31 ) ~ ± w± - f = l^^'^l^) ' 

(ii) Let jV± be the set appearing in Lemma lSTb) (so that 77 is real in the following) . 
Then 

on cA± and for all N G N and a G NQ d+1 ; where the O -symbols are uniform on 
compact subsets of JV±. 

Proof. We drop all ±-indices, set A := dZ — Q*u, and use 

(3.32) d t Z = ( E I V 5 a(Q, H) ) - a(Q, E) = ( E \ d t Q ) - a(Q, E) 

to obtain 

A = {Z' y -{E\Q' y ))dy + {Z^-{E\Q' v ))dy+{Z' v -{E\Q' n ))dr ] 
+ (Z^-(E\QL ! ))drJ, 
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where ( E | Q' y ) dy abbreviates ( S | d yi Q ) dyi+- ■ •+( E | d Vd Q ) dyd, etc. Now, let 
k be one of the variables yj,Jjj, VjiVji J = 1, . . . , cZ, and set \ x := Z' H — ( E | Q' x ). 
Using successively (I3.32p . the Hamiltonian equations (13. 6p . and the almost 
analyticity of a we find 

d t \ K = d«((E\d t Q) - a(Q,E)) - (d t E\d x Q) - (EldtdvQ) 

= ( d x E | V ? a(Q, S) ) - d x (a(Q, £)) + ( V x a(Q, S) | «9 X Q ) 

= - ( VMQ, E)\d=Q)-{ V f a(Q, S) | &*E) . 

Thus, by virtue of (13. 3p . d t \ x = G K , where 

(3-33) d(t,ny,5y,ftri,5ri)G>< = ® ( I ^(Q, | W ) . 

The initial condition (Z, Q, S) | £=0 = ((y \ r) ), y, rj) implies 

Z' y \t=o = rj , Z v \t=o = V ' Zjf\t=o — %rf\t=o = 5 

Q' y \t=o — 1 j Q v \t=o — Qj?-|i=0 = Q77 

and we conclude from diAj, = G K that 



H t ,y,v)-yd , n = y2 G*{r,y,v)dr, {t,y,v) ex>(«) , 

which together with (13.331) yields (i). Part (ii) now follows from Lemma l3~6l □ 
Since we have 

(3.34) Z± o K± = ij} ± , Q ± o K± = x , E ± o K± = g ± , 
on ^ ± we arrive at the following result: 

Proposition 3.9. Let ip± be defined by (I3.25p . Then 

(3.35) <9(W^) {d*l>± + a±(x, V^ ± )) = C(Tj) 
(3-36) ^n.^) (V»V± = O(r^) 
(3-37) flfobua.,ib,)(V^± - **) = 0(r£) 
(3-38) ^.sh^VWi = 0(T£) 
(3-39) ^,^)V#± = O(r^) 

on ^/(4 ; /or JVG N and every multi-index a e Ng d+1 . All O -symbols are uniform 
on compact subsets of j¥±. 

Proof. Lemma I3.8( ii) and T = T o K imply 

d^ x ^ 1l) K* ± (dZ ± -Q* ± u ± -ydr ] ) = 0(T») 
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on K~ 1 (J/±) = JV±. On the other hand f!3.34p shows that, on 

K* ± (dZ± - Q*±LO± -ydrj) = d(K* ± Z±) - (6± o K±)*u± - fc* dr] 

= dvp± — g ± dx + a± (x, g )dt — k^ 1 dr\ . 

□ 

The last corollary of this section summarizes the properties of ip± on the real 
domain, where the weight T± can actually be replaced by S£0±, so that we 
arrive at the desired solution of the problem (I3.1l) fe (l3.2l) . 

Corollary 3.10. (i) There is some real neighborhood, , of <§' ± in Rq x R 2d 
such that, for all (t, x, rj) G ^C R , 

(3.40) $ty±(t,x,7i) > ^yl^^x,^)! 2 , 

(3.41) 5fy + (t,x,7i) = (t,x,r])e£ + , 

(3.42) $h/)-(t,x,ri) = <=> t = 0,x,r]eR d . 
Consequently, 



± 

''R ' 



(3.43) ^ — ^ r± on Jt x 

so that (13.351) - (13.391) hold true on with the right hand sides replaced by 
Oiy((^sip±) N ) . In particular, 

(3.44) dfa,^) (M ± + a ± (x, V^±)) = 0((^ ± f) on Jt± . 
(All O -symbols are uniform on compact subsets of -M^r ■) 

(ii) For all (t, x, 0) G @ + and (5 G N^ +1 ; 

(3.45) d^ x) i> + (t,x,0) = 0, V v il> + (t,x,Q) = X{-t,x,V<p{x)). 
Proof, (i): On account of < K27\f and ATM]) . 

(3.46) ^±-(3x1^) ^ -C(l) ^(x,^)! 3 on Jf±. 

We recall that <A± is a neighborhood of tS^ in IRq" x C d x R d and, hence, 
M±)r := ^± (R^ x R 2d ) is a neighborhood of g± in [Rq x R m . Now, let 
(to, x , rjo) G S ± and let if C (^±)r be a compact neighborhood of (t ,x ,r] ) 
in (~4±)[r. By choosing e > sufficiently small we can ensure that (t,x + 
eQg ± (t,x,r]),r]) G K' C (=^±)ir, for every (t,x,rj) G if and |e| < £ , where if' 
is compact, too. According to (13.461) there exist C,C G (0, oo) such that, for 
all (t, x, rj) G if, 

$hf>±(t,x - e^g ± (t,x,r]),r]) ^ -C \^g ± (t, x - e Qg ± {t, x, rj), rj) | 3 
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Taylor expanding the left hand side of the previous estimate with respect to x 
using (I3.36P and ( I3.38P we obtain 



^ip±>e\Qg ± \ 2 -C"(\Qg ± \ 3 + e 2 \Qg ± \ 2 )-eC No \Qil; ± \ No on K, 

for some N G N, N ^ 2, and C", Cjv G (0, oo). Now, we choose e G (0, 2^77), 
s < e , such that eCV ISt/^ " 1 < 1/2 on if. Then 3t/>± + (1/2)|9^±| > 
(e/2)|% ± | 2 - C"|% ± | 3 on if. Next, we recall from (13T2Tj) and that 
Sg" 11 = on (l' ± . Therefore, we may further ensure that C"\^sg\ < e/4 on K 
by possibly restricting the compact neighborhood K of (to, ^o, Vo) suitably and 
we obtain (I3.40p . Finally, the set is defined as the union of all sets K 
obtained as above for every (t , x , rj ) G S ± . 

Next, we prove (EOIj) and (EQ2]) . First, let (t,x,rj) G S ± . At t = we 
have $^±(0, x, 77) = 3( x \ r\ ) = 0. If (t,x,0) G then we know that 

fc(r, x, 0), 0) G -K"o x {0}, for all r G [0, t], whence ^+(/<r, x, 0), 0)) = 0, 
r G [0,£], due to (13 .4p . Recalling the definition (I3.25P of ip + we see that 
ip + (t,x,0) = 0, which also proves the first assertion of (ii). 

Conversely, assume that (t,x, 7/) G with t > and Qip±(t, x, rj) = 0. 
Then (13T4DD implies that Qg ± (t,x,r ] ) = and (I3T25D and <KWi show that 
(t/, 77) := (fc ± (t, 77) , 77) and (x, 5' ± (t, x, rf)) are connected by a purely real in- 
tegral curve of J^f a± - Then < 3a±(y,r]) < implies $sip±(t, x, 77) > on account 
of (I2.17p . (12.25ft . and (13.27ft . In the minus-case we thus get a contradiction to 
(12.251) showing that there is no (t, x, rf) G with t > and ^sip_(t, x, 77) = 0. 
In the plus-case it follows that Qa + (y,r]) = 0, that is, y G K and 77 = by 
fl2HE]) . Lemma O implies that (t,y,0) G @+, thus (t,x,0) G 

Finally, Part (ii) follows from (j^TSj) . (gJSD , (I3^T|) . and (gSSD -EZ D - □ 

4. The transport equations 

4.1. Formal ansatz for a parametrix. In order to construct a parametrix 
for the conjugated Dirac operator D h y^ we split, roughly speaking, D h y^ 
micro-locally into a plus and a minus part by means of the projections intro- 
duced in (I1.12p . For each of these parts, again roughly speaking, we construct 
parametrices for the corresponding "heat equations" (backwards in time in the 
minus case) and integrate the latter with respect to the time variable. The 
parametrices for the heat equations are obtained as Fourier integral operators 
with complex- valued phase functions. More precisely, our ansatz for the Green 
kernel reads 

fl e{+ ,_ } JO J ® d u=0 [Z7Tn) n 

+ q(x,x- y) . 
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The symbol q additionally appearing here accounts for the elliptic part of Dh,v,ip 
and is constructed in Section below. To find equations determining ip± and 
B± we calculate formally 

oo 

e -#± A (±hdt + a-(-ihV + + a + V) e^ ±/h h " B l 
= ( ± i d t ip± + a ■ (V^± + iVtp) + a + V) B° ± 

oo 

(4.1) + ( ± id t i)± + a • (V(bV± + iVy>) + a + 5±} = . 

In the sequel we fix a smooth cut-off function, x C^°(C 2d ), such that x = 1 
on some small real neighborhood of K x {0}, ^ % ^ 1 on R M , and such that 
supp(x) is contained in some small complex neighborhood of Kq x {0}. We 
assume that x is an almost analytic extension of xtiR 2d ; so that 

\ d (ky&y,Slr,,Qr,)V(y,ri)X{y,V)\ < C N>a My, V) f , (l/.l) ^ M 

for all N G N, a G N^, and suitable constants C^^a G (0, oo). 

Let us suppose for the moment that the matrix- valued amplitudes B± satisfy 

(T ) : df t {B° ± (t, x, rj) - A ± (x, V^±(*, x, rf) + iS?<p(x)) B ± {t, x, rf)} 

= 0(r£), 
B^(0,x,r}) = x(x,T])A ± (x,T] + iV(p(x)). 

From (14.11) we further obtain the transport equations 

(T,W : flf^a.,,) { (Tft + «a • V) i^ 1 

- (±i$V± + a • (V^± + *V<^) + a + V) B v ± } = O n (T^) , 

-D + |t=0 — ~ -D^|t=0 ■ 

If (To) is fulfilled, then the matrix in front of B± in 04. ip can be replaced by 
one of its eigenvalues and we find the eikonal equations 

±id t ^ ± ± y/l + (V x ^± + tV<p) 2 + V = O(rJ) , 

which, according to the definition f !2.16p . are equivalent to the problems (13.11) 
solved in Section [3J Again the error terms 0(T±) in the transport equations 
cannot be avoided because the transport equations are complex-valued. They 
do, however, not destroy the WKB construction as we shall see later on in 
Proposition 14.81 
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4.2. Solution of the transport equations. In the rest of Section |4] we as- 
sume that V fulfills Hypothesis II. 1\ x* and fulfill Hypothesis ll.2[ and that 
ip and ip± are the functions given by Propositions 12.11 and 12.21 and by (I3.25p . 
respectively. All (9-symbols are uniform on compact subsets and the variable 
r\ will be real. 

The transport equations (To), (Ti), . . . are solved by means of a strategy we 
learned from [llj. Following this strategy we have, however, to keep track of the 
error terms and to put some factors % and some additional minus signs in the 
right places. For we consider a "heat equation" with different time directions 
for the plus and minus parts of the symbol rather than the usual Dirac equation 
whose scattering theory is discussed in [TT] . For the convenience of the reader it 
thus makes sense to give a self-contained discussion of the transport equations. 
As a first step we introduce gamma-matrices 

7o := Qt , 7j := -a a>j , j = l,...,d, 

so that 

(4.2) (7o) 2 = l, (7;) 2 = -l, j = l,...,d, 

(4.3) { 7m)7i/ } = 0, 0^/2<u^d. 

Furthermore, we set 

:= ±id t , df := -d x . , 
II± := ± yfi + (V X V± + *V</?) 2 , Uf := d Xj ij ± + id x . <p , 

n± : = (n± ...,n±), 

where j = 1, . . . , d, and 

d d 
^:=E^^' n ± :=^ 7 ,Hj on ^ ± . 

We recall that the sets have been introduced in (13.290 . We transform the 
transport equations into a new sequence of equations on jV± given as 

(K , \d(tM^,v) (( ft± - 1) Bl + i~d±B^) = O n (T») , N E N , 
\^,^)( ft± + 1) ^ 5 ± = 0»FZ) , N e N , 
where z/ G N and S^ 1 := 0. 

Lemma 4.1. If B v ±, v G No, satisfy the first equation in (Ki), (K2), . . . on 

some neighborhood, jY± C ^44, 0/ <£± m Rj x C 1 x R d ; i/ien i/iey satisfy the 
transport equations (Ti), (T 2 ), ... on too. 
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Proof. Let v G N. Multiplying the first equation in (Kj,) with «o and using 
«q = 1 we obtain 

d^ x , v) {(Td t + » • V) i^" 1 + (itf + ^ 

(4.4) - a • (V^± + iVp) - a - ^) 5 ±} = 0jv(r£) . 

In view of the Hamilton- Jacobi equation 52 ^ 9 ^ (±id t i^± + Tl^ + V) — 
O n (T±), Af 6 N, we see that (JO]) is equivalent to (T„). □ 

In what follows we abbreviate 

(4.5) <v := % - [ T „ , lv ] = z T „ lv , (rot ± 11% := (c£n±) - ) , 

for all //, f = 0, . . . , d, fi ^ u. We introduce the following matrix-valued partial 
differential operator on <A±, 



(4.6) L± := ±i{d t ,U±} + ^2{d Xj ,Uf}-i ^ o> (rot ± n d 
On account of (j4.2j) . ( I4.3p . and ( 14.5)1 we can re- write L 1 * 1 as 

(4.7) l± = i d t > nj } + £ 7, 7, (<t o n± + nj a±) = , n*} . 

Lemma 4.2. Let jV± be the neighborhood of S± in Rq x C d x R d introduced 
in (13.291) and suppose that the matrix-valued function R G C°°( e A±, Jzf (C d *)) 
sate/ies 

(4-8) flft^sh,,)^ = 0(T£) , iV e N , a G 

Furthermore, let c G C°°(C d x R d ) satisfy 
(4-9) ^.v.;,,,,^ = 0(19*1"), iV G N , a G 

T/ien i/iere exzsi 5 ± G C°°(^4, j£? (C d *)) suc/t tfiai 

(4-10) flK^-^fl* -R) = 0(lf), dejB% =0 - c) = 0(\Qx\ N ), 



(4.H) % H *, a . f ,)^S ± = 0(I±), 

/or ail iV 6 N, a 6 N^ +1 ; and /3 G N|j d . G C°°(^, J^(C d *)) « another 

solution of (TCTUj) and (Oil) , toen <9 ( Q t j^^., - C ± ) = 0(r£), /or a// 

iV G N and a G N^ +1 . 
Proof. By definition we have 
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where 

m ± -.= ±i(d t n±) + div :c n ± - i <v( rot±ir V- 

We further set 

W ± := -(±2in^)~ 1 M ± , S ± := (±2in^)~ 1 R , 
and consider the maximal solutions, B^, of the ordinary differential equations 

(4.12) dtB^ViV) = W ± (t,Q ± (t,y,r ] ),r ] ) B ± (t,y,r ] ) + S ± (t,Q ± (t,y,r ] ),r ] ) 

defined for (t,y,r)) G *A± with boundary condition B(0,y,rj) = c(y,T)). We 
further introduce the vector field 

Z± := V 5 a±(x,5f ± ) • V x 

on < ^ ± . Then an appropriate restriction of : ©(k 1 * 1 ) — > C d is equal to the flow 
of the real vector field Z± := Z±+Z±, because g ± (t, Q (t, y, rj),r]) = H ± (t, y, 77), 
for all (t,y,rj) G ©(k 1 * 1 ) such that (t,Q ± (t,y,r]),T)) G ^ >± . Setting 

B ± (t,x,r ] ) := B ± (t,k ± (t,x,r ] ),r ] ) 

we thus have 

(d t B ± + Z ± B ± )(t,Q ± (t,y,r ] ),r ] ) = j t B ± (t,Q ± (t,y,r ] ),r ] ) 

= W ± (t, Q ± (t, y, 77), 77) B ± (t, y, 77) + S ± (t, Q ± (t, y, 77), rj) , 

because k ± (t,Q ± (t,y,r)),r]) = y. In view of (±2«ng) -1 2II ± = V£a±(x,V x ip±) 
it follows that 

(±2zn±)- 1 (L ± B ± -R) 

= (d t + V e a±(x, VsVi) • V x - W±) 5 ± - S± 

(4.13) = -Vsa ± (x,g±) ■ V^ ± + (V € a ± (x, V^±) ~ V f a ± (x, • V,5 ± . 

From now on we drop all sub- and superscripts ± in the existence part of this 
proof since they do not play any role anymore. By virtue of (13.361) we know 
that ^,^,3,,^^ ~g) = 0(T N ), whence 

(4-14) ^U„„„ /: (V € a(s, V x ^) - V 5 a(a; ) g)) = 0(T N ) , 

for all iV G N and a G Nq^ 1 . Furthermore, it is clear that all partial derivatives 
of B are bounded on compact subsets of jY . To study d^B we write 

d^B = (d y B)(t, k, 77) d^k + (d ¥ B)(t, k, 77) d^k . 
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Here we know from Lemma I3T71 that d^ t ^ ^ v ^d x k = 0(T N ) and it suffices 

to find a similar bound on dyB. To this end we differentiate the differential 
equation (14. 12|) to get 

d t d v B = W(t, Q, rj) d w B + (d x W(t, Q, rj) dyQ) B + (d^W (t , Q, rj) dyQ) B 
+ d x S (t, Q, rj) dyQ + d x S (t, Q, rj) dyQ , 

which on account of fl3~T5]) . (13381) . (Oj) . and \<3Q\ 2 ^ 0(1)?, shows that 
(4-15) &( my ^ v) (d t d v B - W(t,Q,rj) d v B) = 0(f N ) . 

Since B(0,y,r]) = B(0,y,r)) we also have d^ %fl) d ¥ .B| t= o = d^^dyC = 
0(\^sy\ N ). Therefore, we first obtain from Duhamel's formula and T(s,y,r]) ^ 
O(l)f(t,y,r}),se[0,t],tb&t 

(4.16) d w B = U t , O(\$y\ N ) + [ U t , s O(f(s,y, V ) N )ds = 0(f(t,y, V ) N ), 

Jo 

where U tjS fulfills d t U ttS = W(t,Q(t,y,r)),T])U tjS , U SiS = 1. Now, suppose we 
have shown that dL yQyv} d w B = 0(f N ) is valid, for all /3 G Njj d with < n G 
N and let a G Nq^, |a| = n + 1. Then we obtain 

(4.17) d t df m>Qy>v) d v B - W(t, Q, V ) 9fc,*y„)<hB = 0(f N ) 

from (14.151) . and we again conclude that d^L^y^dyB = 0(T N ). Using the 
differential equation (14.151) and the usual bootstrap argument to include the 
time derivatives we further see that d" t ^ y ^ y ^dyB = 0(T N ) , for all N G N and 

a G No rf+1 . Altogether we arrive at (I41ip and from (I4.13P and (14. 14[) we infer 
that (I4.10p holds true also. 

Now, suppose that C 1 * 1 is another solution of (I4.10p and (14. lip . Then we have 
{±2tU^)- 1 L ± {B ± - C ± ) = G ± , where d^^^G* = 0(T»), or (compare 

[d t + Z±- W±) (B ± - C±) = G ± 
+ (V c a ± (x,^) • V* - (V € a ± (x, V^±) - V € a ± (x, g*)) ■ V X }(B ± - C ± ) . 

Let denote the right hand side of the previous identity. Setting E ± : = 
(B ± - C ± )(t, Q ± , rj) we then have 

d t E ± -W ± (t 7 Q ± ,r 1 )E ± = 17), ^ ytrl) E% =0 = 0(\%y\ N ) , 

where duvtxSSxrj)^ = 0(T±). By the same induction argument as the one 
used above to discuss dyB we infer that d? t ^ y ^ y ■> E ± = 0(T±), which implies 
dft^jB* ~ C±) = 0{T» ), for N G N, a G N** 1 . □ 
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The next lemma will be applied to the choices f ± = id ± B± and u ± = 

Lemma 4.3. Let the matrix-valued functions f ± ,u ± G C°°{J/±, (C d *)) sat- 
isfy the equations 



(4-18) fi&i te ,a..,)(n ± + i)/ ± = e>(r£), 

(4.19) 9 ( V^ ) (L ± n ± + 9 ± / ± ) = 0(r£), 

(4-20) flS.ito.ar,,)^ = C(r±), 

(4.2i) ^;...,,,, ;: ^.r = o(r£), 



/or every JVe N and a G Nq , and 

(4.22) d^ x>v) ((n±(0, x, ry) - 1) ^(0, x, 77) + f±(0, x, 77)) = 0(\$*x\ N ) , 

/or a// (0,37,77) G iV G N, and (3 G Nj] d . 27jen /u//i//s £/ie following 
equations on JV±, for every N G N and a G No d+1 , 

(4-23) (( ft± - 1) « ± + / ± ) = 0(1^) , 

(4.24) ^^(n* + 1) ffM = o(r£) . 

Proo/. On account of (TOD and (TOD we have (fl ± ) 2 = (Il±) 2 -(V x ip±+iV(p) 2 = 
1, which together with ( 14. 7p gives (IT* 1 — 1) L ± = L ± (IT — 1). From this we 
infer that 

l ± ((ri* - 1) m ± + /±) = (n ± - 1) l ± u± + (^n* + fi ± 9 ± ) f± 

= (fi ± -i)(L ± u ± + d ± f ± ) + d ± (u ± + i)f ± , 

thus d^ ^Z^ ((fl ± -l)M ± +/ ± ) = £>jv( r ±)- Together with Lemma Oand 
(T4221 this shows that (14T23]) holds true. The identity J4T24) follows from (TCT9]) 
and (1423]) because (n ± + l)9 ± M ± = L± m ± + 9 ± / ± -9* ((n* - 1) m ± + / ± ) . □ 

By virtue of Lemma we may now define P4 G C°°(„4±, j£?(C d *)), f G N , 
by successively solving the differential equations 

(4-25) (L* B± + id^Bl' 1 )) = 0{lf) , := , 

on o/Kj- with initial conditions 

(4.26) B° ± (0,x, V ) := X (ar, r?) A ± (n ± (0, x, 77)) = v) ^(v + W<p(x)) , 
for (0, x, 77) G J/±, and 

(4.27) ^| t=0 := -(2n±)- 1 a (79 + J B;- 1 + z,9- J Br 1 )| t=0 , ^N. 

Here x is the cut-off function introduced below (14. ip . We summarize the results 
of the previous constructions in the following proposition. 



28 



Proposition 4.4. The matrix-valued amplitudes B± G C°°(o/4±, Jzf (C d *)) c?e- 
/med fry fl4.25p - fl4.27l) so/ve the original transport equations (To), (Ti), (T 2 ), . . . 

on jV± and 

(4-28) flg,^,,)^ + 1) ~d ± Bl = 0(T!*) , N G N , a G N^ 1 . 
Moreover, 

(4.29) S(V,^)^ S ± = ( r ±) . ^ e N , a G N** 1 . 

JTie supports of B±, v G No, are compact and contained in some fixed compact 
neighborhood of 

Proof. We argue by induction successively applying Lemma [4.31 For v = 0, we 
set f ± = so that flggp is satisfied trivially and ( KWf is just ffOo) . The 
initial condition f)4.22p is implied by f)4.26p . By Lemma [4.31 we see that (I4.28P 
is valid, for v = 0, and that (n± - 1) S£ = 0(r£), which implies 

5 ("5Rx,^) AT (Vx^± + iV^)5° = 0(rf), which is (To). Next, assume that 
n G N and that B u ± fulfills (T„), (jOgjl . and p5| . for every z/ = 0, . . . ,n. 
Then (Qgp - dOIjl are fulfilled with / ± = i^SJ and m ± = by (Qgjl 

and the definition of At t = we have 

(n ± - 1) ^ +1 | t=0 = (-2n 5 t )- 1 (fi± - i) a (z9 + i?: + z5-5!!)| t=0 

(4.30) = (-2n 5 t )- 1 (n* + a ■ n* - a ) (*<9 + ^ + id~B n _) | t=0 . 

Furthermore, we observe that (I4.28P with u = n yields 
(4-31) dfL^W (<^)U - (a H± - « ) (<^)U} = 0(19*1") , 
because T(0, x, 77) = \5s(x, rj)\ 2 — ( 1 9ft?7 ) +9 : (7/ 1 x ) = |Sx| 2 , for real 77, where 
(4.32) n ± | t=0 = 77 + iV^(x) 

is independent of the choice of the ±-signs. Combining fl4.30p - fl4.32p we arrive 
at 

= ^^{(-mtr 1 ((nf + n+^d+Bi) + (lit + nzXiSrBl)) \ t=0 } 

In the last step we used 

(4-33) n + | 4=0 = Vl + (r/ + *V^) 2 = -IL7| t=0 . 

In conclusion we see that (I4.22p is satisfied, too. Again by Lemma 14.31 we 
deduce that ( jOgJ) is fulfilled, for v = n+1, and that ^ ^ ^ ^((^±-1) Bl +1 + 

id ± Bl) = 0(T%). By virtue of Lemma IP we conclude that Bl +1 fulfills the 
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differential equation in (T n+ i). The initial condition in (T n+ i) is satisfied 
because of ffl~2Tj) and Q£ 



Finally, we recall that the estimates ( I4.29P follow from Lemma 14.21 The last 
assertion on the supports of B± is clear from the constructions in Lemma I4.2[ 
the fact that the values of t are bounded on <A±, and the fact that the initial 
conditions B±\ t =o are supported in supp(x)- □ 

In order to calculate the leading asymptotics of the Green kernel of D^y at 
(x*,y*) we have to compute the value of -B°(r, x, 0), for (r,x,0) G S> + ', more 
explicitly. We recall that, by definition, (r, x, 0) G ^ + implies that there is 
some y G K such that y, 0) G K , t G [0, r], and Q + (r, y, 0) = x. 

Lemma 4.5. Le£ 5° fre a solution of (T ) as m (I4.25P and (I4.26P . Lei 

(r,x,0) G J> + and Zei U(-,y) : [0,r] -> Jzf(C d *) be the solution of 

d TT , x «a W(Q+(t,y,0)) rr , v rr/n , 

on [0, r], where Q + (r,y,0) = x. Then 

(4.34) 5j(r, *, 0) = (-ny))^ x (y 0) ^ A+ 

(-V(x))V2 de t[ci J/ Q+(r,n,0) 11 - 



Proof. On account of Proposition ^. 9l we have x, 0) = 0, V x ip+(t, x, 0) = 0, 
thus n+(t,x,0) = y/l - V(p(x) 2 = -V(x) and n + (t,x,0) = iVp(a;), for all 
(t,x,0) G We further deduce that d t U.^(t,x,0) = 0, drv x ll + (t,x,0) = 

iAp(a;), <9+n+(t,x,0) = -id Xi d Xj (p(x), d+U+{t,x,0) = 0, dfli+(t,x,0) = 
d x .V(x), where i,j G {1, . . . , off, thus 

-2 ^ a ^ ( rot+ n+ W*> ^, 0) = 7o 7 ■ (- W(a;)) = a ■ VV{x) , 

for all (i, x, 0) e @ + , where 7 := (71, . . . , j d ). For (i, y, 0) G the differential 
equation (14.121) determining , for the choice 5* + = 0, thus reads 

(4.35) d t B + (t,y,0) = 2W{Q+{t, y, 0)) °)> 

for £ G [0, r], and we have the initial condition -B+(0, y, 0) = x{Vi 0) A + (2Vy2(y)). 
Using the ansatz B^_(t,y, 0) = b(t,y) U(t,y) Bl(0,y,0) with a scalar 6 it thus 
remains to solve 

y) = ^(g+^^oj) ^^ y) ' * G [0, r] ' ^ (0, y) = 0) ' 
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Using Liouville's formula for the equation d t Q + (t, y, 0) = F(Q + (t, y, 0)), where 
F(x) := V p H(x, Vif(x)) = -V- 1 ^) Vy?(x), x G K , that is, 

d t det[d y Q + (t,y,0)} = divF(Q + (t,y,0)) det[d y Q + (t, y, 0)] , (t,y,0) G S> + , 
where divF = — V^ 1 A(p — V~ l W ■ F, it is, however, elementary to verify that 

d t {(-V(Q+(t, y, 0)))" 1 / 2 det[d y Q + (t, y, 0)]" 1 / 2 } 
A<p(Q+(t,y,0)) 



2V(Q+(t,y,0)) 



{ (-V(Q + (t, y, 0)))- 1 / 2 det[d y Q + (t, y, O)]- 1 ^} . 



We deduce that b(t,y) is equal to the term in the curly brackets {• ■ ■ } times 
(-V(g+(0,?/,0)) 1 / 2 = {-V{y)) x ' 2 (so that b(0,y) = 1) and the formula 
follows. □ 

In the one-dimensional case the formula (I4.34p for the solution of the first 
transport equation can still be written a bit more explicitly. 

Lemma 4.6. Assume that d = 1, let B^ be a solution of (T ) as in (I4.25P and 
(Q5]l . and let (r, x, 0) G # + . Then 

B °(rx0)= t- v (y)) 1/2 x(y,Q) 

+[ ' ' J (-V(x))VHQ+y y (T,y,0)V* 
(4.36) ■ ( cos(0(t)) 1 - i sin(tf(r)) on) A + (itp'(y)), 

where Q + (r, y, 0) = x and 

nt) := / 9T/ , n+ , T7rzds, te[0,T\. 

Jo 2V(Q + (s,y,0)) 

Proof. On account of Lemma [4.51 we just have to compute the solution of 

which is given by U(t, y) = cos($(i)) 1 — i sin(i?(t)) oti, t G [0, r]. □ 

Remark 4.7. We know that n + (r, x*,0) = — ct ex ■ 2V<£>(x*) + cto (— V (#*)), 
where V <£>(£*) = w(r). Then the identity (Il + — 1) B9(t, x+, 0) = implies 

B+(t,x^0) = (i + ^n + (r,x„0)) B+(t,x*,0) 

= (-\/(^))A( ra (r))a 5 +(r,^,0). 

In the following we verify directly that taking the hermitian conjugate of 
M(a;*,j/*) := (-V r (x Jr ))A + (^(r)) a U(t) (-V( yir ))A+(iw(0)) gives the same 
expression as interchanging the roles of x± and y*. 
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Notice that the Hamiltonian trajectory in {H = 0} whose position space 
projection runs from X* to y* is given by (~) (i) := (_ 7 ) (r — t), t G [0, t). Thus 



M(y*,x*) = (~V(y,))A + (-tw(0))a U(r)(-V(x^))A + (-tw(r)) 
solves 

?7(t) = tot ■ w(j(t)) U(t) , t G [0, r] , w := -(2K) _1 W . 



where U solves 
d 
dt 

Hence, in order to verify M(y+,x+) = M(x±,y±)* it suffices to show that 
a U(r) = U(t)* a . To this end we recall that U(t) is given by the Dyson 
series U(r) = X^ =0 4(t), where 



I n (r) := / ia-w("f(t n ))---ia-w(j(t 1 ))dt 1 ...dt n , 
JtA„ 

t A n := {0 ^ ti ^ . . . ^ t n ^ r} denoting the n-dimensional standard simplex 
scaled by r. On the one hand, {«o, %} = and a* = ctj, j — 1, . . . , d, implies 



4(t)* a = a ict- w(^(ti)) ■ - ict- w(^(t n )) dti . . . dt n . 
JtA„ 

On the other hand the substitution s x = r — t n , . . . , s n — r — 1\ turns the latter 
expression into 



a ict- u>(7(s n )) ■ ■ ■ ia ■ 10(7(51)) ds x . . . ds n =: a I n (r) , 

J T A n 

where £/(r) = £~ 7 n (r). O 

4.3. Approximate solution of (±hd t + Dh,v,ip) u ± = 0- Finally, we put 
the results of Section [3] and the present section together. To this end we pick 
smooth cut-off functions, g± G C°°(IRq x R 2d ), such that g± = 1 in some 
neighborhood of $± and supp(^±) C JV± fl (Rq x R 2d ) and such that all partial 
derivatives of any order of q± are uniformly bounded. Furthermore, we define 
B±{ -;h)e CH^o x R M , ^(C d *)), h e (0, 1], by 

00 

(4.37) B ± (t, x, m h) = J2 hV 9 ( h / £ ») B±(t, x, V ) Bl{t, x, V ) , 

for (t, x, 7]) GKj x R 2d , h G (0, 1], where 9 G C°°(R, [0, 1]) equals 1 on (-00, 1] 
and on [2, 00). (B± is set equal to zero outside its original domain of definition 
JY±.) If e u \ sufficiently fast, we have 
(4.38) 

AT 



d (t,x, v ) [B±(t, x, r);h)-^2 h " 0±(*> x ; ^) 5 ±(^ x ; V) 



sup 

4- v — U 

RjxR 2d 



^ C/VfY h N+1 , 
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for iV G N, a G Ng + , h G (0, h a ^], and suitable constants CV^, ft-a.jv £ (0, oo). 

Proposition 4.8. There exist compactly supported matrix-valued functions, 
r±(-;h) G C°°(IRo x R 2d , Jz? (C d *))> such that 

(±h d t + D htVtV ) {e^' h B±) = f± , lie (0, 1] , 

and, for all N G N and a G NQ d+1 ; t/iere zs some CVq, G (0, oo) swc/j t/iat 

(4.39) sup \\d% iXiV) r±(t,x, m h)\\ ^ C N , a h N , h G (0, 1] . 

(t,x,r;)elR f |xlR 2d 

Proof. Since all terms in (14. ip corresponding to the different powers of h are 
equal to some smooth matrix-valued function on jY± whose partial derivatives 
of any order are equal to 0(T±), N G N, and since T± ^ 0(1) Qip± on the real 
domain, 

oo 

(±hdt + D h ,v iip ) [e^ h B ± ] =J2 h " 0{h/e v )^' h O({^ ± ) N ) 

oo 

+ '$2 h" 9{h/e v ) e^ ±/h [(±d t - iot ■ V)g ± ] B v ± l 

u=l 

oo 

(4.40) +J2 hU (W^-i) " 9(h/e v j) e^ /h (±d t - icx ■ V) [f?± . 

u=l 

Consequently, all partial derivatives of the first term on the right hand side of 
(I4.4UI) are of order 0(h°°) because 

(4.41) e -^±/h ^ ± ) N <: N\h N , N G N . 

Furthermore, $fip± > on supp(p'_j_) and all partial derivatives of R u := [(+d t — 
ia ■ V)q±}B v ± are locally bounded, whence df t>x>r)) [ e i ^ h R v ] = 0(h°°). All 
partial derivatives of the third term on the right hand side of (14.401) are of 
order 0(h°°), too, since - d(h/e v ) = 0, for all h G (0, e„), u G N. □ 

5. A PARAMETRIX FOR D h y^ 

Given fc,m G R, we write 6 G S k ((£) m ; J? (C d *)), for some map 6 : R 2d x(0, fo ] -> 
Jzf(C d *), if > 0, b{-;h) G C°°([R 2d , Jzf (C d *)), h G (0,/i ], and, for all a G N 2d , 
we find h a , C a G (0, oo) such that 

\\d^ b(x^;h)\\ ^ C a (O m h-\ x,£GR d , he(0,h a ]. 

We further set 

S-°°((£}-°° ;J Sf(C d *)) : = p| S- fe ((0^';^(C d *)). 

fcGN 
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In what follows we work with the semi-classical standard quantization of matrix- 
valued symbols b G S k ({£,) m ; J2?(C d *)) determined by the oscillatory integrals 

Op h (b) f(x) := Je^ x ^ h b(x,Of(y)-^- d , /G^(R d ,C d *). 

Let k 1 ,k 2 ,m l ,m 2 G R. We recall that, for b G S kl ~^(C d *)) and c G 

5 fc2 ((0 m2 ;=^(C d *)), the symbol, b# h c G S* 1+ta (((> mi+ras ; J^(C d *)), of Op h {b)o 
Op h (c) has the following asymptotic expansion in S kl+k2 ((^) mi+m2 ; J?f (C d *)), 

&# fcC (x, £;/») = &(*, v) c(y , 1 — - Z) 4n ^ 6 ) ^ ^ ^ ^ ' 



If the symbol Dv >v (x,£) G S°{(0] ^{C d ')) defined in f[2TToD were invertible, 
for every G R M , we had a well-known asymptotic expansion, q(x,£) x 

Yl^Lo h u Qv(x,Q, of the matrix- valued symbol of the inverse operator. We can, 
however, write down this asymptotic expansion formally and determine q u (x, £) 
at every point (x,£) G R 2d where Dy ttp (x, £) is invertible. We proceed in this 
way and pick some cut-off function x E C^°(R 2d , [0, 1]) such that x = 1 in a 
small neighborhood of K x {0} and supp(x) C {x = 1}, where % has been 
introduced below Then it turns out that g v (1-%) G S°{{g)~ v ~ l \ £?{C d *)). 

Let g be a Borel resummation of Yl'^Lo ^ 1" 0- ~~ X)> so that 

N-l 

? - 5>" *> t 1 ~ x) G ^((O"^" 1 ; ^?(c*)) , at g n . 

i/=0 

Then D h y jlf o Op ft (g) is a pseudo-differential operator whose symbol has the 
asymptotic expansion 

oo 

D v ^# h q{x,£) x 1 - x(x,0 + ^2h v r v (x,^) . 

Here each error term r v contains some partial derivative of % whence supp(rj,) C 
{x = 1}, for every v G N . Setting 

q := g#fc(l - %) 

we thus have 

(5.1) D v ^ h q- (l-xjenr 00 ;*)). 

Next, we define an operator V h : y{R d , C d *) -»■ •*"(R <i , C 1 '), 

dydrjdt 
(2nh) d h 



(V h f)(x):= T M / e^^^-^l^/^tt^x,^/.)/^) 
(5.2) + Op h ( ff ) /(x) , x G R d , / G ^(R d , C d *) . 
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Here the integrals in the first line are effectively evaluated over some compact 
set so that Vh is obviously well-defined. Furthermore, it is clear that Vh can 
be represented as an integral operator with kernel 

V h (x,y)= tt/X 6 ^'^''"^^^'^^ 5 ^^^ 



fle{+,-} 

(5.3) + q(x,x-y) . 



We recall that q(x,y — x), with q(x,y) = (J r h ~ 1 )^ y q(x,y), is the distribu- 
tion kernel of Op h (q). Here we normalize the (component-wise) semi-classical 
Fourier transform as 

f( v ):=(F h f)( V ):= [ e-^/ h f(y)dy, V G R d , f G L\R d , Y) , 
where V is C d * or J2?(C d *). Integrating by parts by means of the operators 



1-ih (Vt/0± -y)-V v 
l + |V^ ± -y| 2 

and using the fact that B± is compactly supported it is easy to see that 
\\d^V h (x,y)\\ < Cwh-W-W (x)~ N (y)- N , x,y G R d , 

for all N G N, a, (5 G Nq, and suitable constants CN, a ,p G (0, oo). 
Theorem 5.1. (%) There is some f G S^ 00 ^) -00 ; Jz^(C d *)), such that 

D h , v , v V h = l-Op h (f). 

(wj T/iere exist h G (0,1] and smooth kernels, K h G C°°(R M , if (C d *)) ; h G 
(0, /io], snc/i that 

(5.4) ||^7^(z,y)|| < C7 w /i JV (x>- JV (y>- JV , x,y E R d , h G (0, /i ] , 
/or a// iVGN, «,/?G Nq, and suitable constants Cn,o.,p G (0, oo), and snc/i that 
D? v (x, y) = e-^W-^W)/" y) + 7^(z, y)) , x ^ y . 
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Proof, (i): By Proposition l4~8l we have, for / G C£°{R d ,C d *), 

dr] dt 



POO f 

D h,v, v Vf e^*M h B t {t,x, m h) fa) 



m+,-} 



(2nh) d h 



p poo 

- E / / dt{& { *» ),h Bfax^h)) fa) 

u^ r , i JR d JO 



tte{+,-} 



poo p 



»e{+,-} 



dr] dt 
(2nh) d h 



(it C?7/ 

(2vr/i) d 



c??7 



(2tt/i) £ 



(5.5) - Op fc (fi)/(x) 



where A + + A~ = 1. We recall that the integral appearing here in the first line 
is effectively an integral over some compact set. In particular, the boundary at 
oo does not give any contribution to the integral / °°9t(- ■ ■) dt. Moreover, we 
put 

poo J± 

fx(x, m h) U e-^/ h f t (t,x, V ;h)j. 

In view of (I4.39P and since f has a compact support it is clear that fj is a 
symbol in S' _00 ((?7) _00 ; Jz? (C d *)). In fact, if some derivative d? xr} \ is applied to 

fx the inverse powers of h obtained by differentiating the phase e~ % ^ x 'l h are 
compensated for by derivatives of fj, which are of order 0(h°°). On account of 
(15. ip we further have 

DH,v, v Op h (q) f(x) = [ e^^il-xix^^f^J^-OM^fix), 

for some f 2 G S 1-00 ^) -00 , Jzf (C d *)), which proves (i) with f := fi + r 2 . 

(ii): For sufficiently small h > 0, ||Op fe (f)|| ^a) ^ 1/2, and (l-Op^f)) -1 = 
Op ft (c), for some c G ^(l; J&f (C d *)). Now, Op h (c) (1 - Op h (f)) = 1 is equiva- 
lent to Op h (c-l) = Op fc (c) Op A (f) which implies c-1 G S^ 00 ^)" 00 ; Jzf(C d *))> 
because f G S^ 00 ^)" 00 ; JSf (C d *))- Of course, V h Op h {c) = V h + V h Op h {c - 
1) and the distribution kernel of Op h (c — 1) - let us call it K-h - fulfills 
\\d^d^JCh(x, y) || ^ CN, a ,p h N (x — y)~ N . Therefore, using the remarks preceding 
the statement of this theorem it is easy to check that IZh '■— Vh Op fo (c — 1) is an 
integral operator with smooth kernel - again denoted by the same symbol - such 
that (J53D holds true. Using D h>V;tp (V h +K h ) = 1, and D hjV)tp = e v/h D hy e~ v/h , 
where ip is bounded with bounded partial derivatives of any order, we conclude 
that e~^ h {V h + K h ) e^ h = D^ v . D 
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6. Calculation of the leading asymptotics 

In this section we calculate the asymptotics of the integral kernel (I5.3P of the 
operator Vh defined in (15. 2 p at the distinguished points x+ and y* fulfilling 
Hypothesis 11.21 On account of Theorem 15.11 this will complete the proofs of 
our main Theorems 11.31 and 11.51 As in the statement of these theorems we let 
'■ [0, t] — > R 2d denote a smooth curve solving (ll.lOp and satisfying (11. lip 
such that 7(0) = y± and 7(7-) = 2* and set 

v» := ^7(0) , := f^r) . 

Proposition 6.1. Let d E N. As h E (0, 1] tends to zero, 

(61) 1 / h_\ ^ (1 + gW) (-V{y*)Y 12 Ufa y*) A + (*Vy^)) ^ 



'-V( Xi ,)y/2 Jdet 



~ v y* 
v Xir id v Q + (r,y±,0) 



where we use the same notation as in f !4.34j) . 

Proof. By Theorem I5.1( ii) it suffices to consider only the kernel Vh- First, a 
standard argument shows that the distribution kernel q(x,x — y) of Op h (g) 
in (15. 3p does not contribute to the asymptotic expansion in (16. ip . In fact, 
(x — y) 2N q(x, x — y) is the inverse Fourier transform of h 2N Af q(x,£) at x — y, 
where A^q(x, £) is absolutely integrable with respect to £, for large iVGN. 
Next, we consider the integral 

At t = we have ^-(0, x±, r/) — (rj\y+) = ( r\ \ x+ — y+ ). Since x± 7^ y+ we can 
thus show by integration by parts with respect to r] that 

nJ4>-{t>x*,v)/h-i{n\y*)/h B ( f r h \ dr l dt _ r)(h°°) 
e u„{t,x*, V ,n) ^ h)dh - ^{n ), 

provided that e > is sufficiently small. Since Qip-(t, x*, r/) > 0, for t > 0, 
by ( EH) , it further follows that f™f Rd e^-/h-i(v\y*)/h B_drjdt = 0(h°°), for 
every fixed e > 0. Consequently, = 0{h°°). 

Finally, we treat the integral 

l + {x*,y*) - J q j^e B + {t,x*, V ,h) 

which is the only term contributing to the asymptotic expansion. We shall 
apply a complex stationary phase expansion with respect to the d + 1 variables 
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(t,rj). The critical points of the phase are given by 

(6.2) = M+fax^r]), 

(6.3) = W n 4>+(t,x±,vi) - y ± . 

To find the asymptotics of I+(x+, y±) it certainly suffices to determine all critical 
points (t,rj) with x*, rj) = 0. We know, however, that $sil) + (t,x+,r)) = 

implies t = or (£,£*, 77) G S> + . As above, we infer by integration by parts 
that J £ J Rd e ^+/h-i(v\y*)/h B+ drjdt = 0(h°°), for some sufficiently small e > 0. 
The only critical point (t, rj) such that (t, x*, rj) G S> + is, however, given by 
(t,rj) = (t, 0). The method of complex stationary phase [6] thus implies that 

(2tt/i)^ ^+( T > x *<°y h B° + (r, Xin 0) (f + 0(h)) 



(6.4) J+(z*,y*) 



Thanks to (13.4=51) we know that t/> + (t, a;*, 0) = = dfip^TjX+jO) and differen- 
tiating the identity Q + (t, fc + (t, £*, 77) , 77) = we obtain 



d y Q + d t k + (t,x*,r]) + dyQ + d t k+(t,x±,r]) = -d t Q~ 



d y Q + d v k + (t,x*,r]) + d ¥ Q + d v k + (t, Xin rj) = -d v Q + , 

where all derivatives of Q + are evaluated at (t, k + (t, x*, 77), 77). Inserting (i, 17) 
(r, 0) we infer that 

(6.5) d y Q + (r,y ic ,0)v yi , = v Xir , 

(6.6) d y Q + (r,y it ,0)\d rj Vr l '^+{T,x ir ,0) = i d ri Q + {r, y+, 0) , 
The determinant appearing in (16. 4p is thus equal to 



(6.7) det 



-v. 



d y Q + (T,y„0)-Hd v Q + (T,y^0) 



Multiplying the determinant (16. 7p and the one appearing in (14.341) . which can 
be written as 

det[^Q + (T,y*,0)] = det(j dyQ+ ^ 0) 

and using x{y*i 0) = 1, (16.51) . and (16.61) we arrive at (16.1 p . □ 
Next, we re- write the formula ( 16. ip in the one-dimensional case. 
Proposition 6.2. For d = 1, we have, as h \ ; 

D _i r v = 1 + Q(*0 

(6.8) • ( cos(#(r)) 1 - i sin(0(r)) e*i) (-V(y*)) 1/2 A(zV^(j/*)) , 
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Proof. In the case d — 1 the formula (16. ip reduces to 

Here w Xi and t>^ have the same sign and |t> 2 | = (1 — V 2 (z)Y^ 2 /\V(z)\, for 
2 = so that the factors |V(.z)|, z = cancel each other in the 

denominator. Moreover, we have already calculated Z7(r, y*) in the proof of 
Lemma 14.61 □ 

In more than one dimension we can evaluate the determinant in (16.11) more 
explicitly. To this end we first prove the following lemma. 

Lemma 6.3. Let X be the position space projection of the Hamiltonian flow 
associated with H as defined in (12.111) . Then the following identity holds, 

(6.9) id,Q + (r, y„ 0) = d p X(r, y+, w(0)) , 

where w(0) = V<^(y*) is the initial momentum of the Hamiltonian trajectory 
from y* to x* as in the statement of Theorem M.SX 

Proof. By Lemma Owe have Q + {t,y*,0) = X(t,y„,0) = j(t), t G [0,r]. We 
set pit) : = (7(t),cj(t)) = (7(f), V<p{l{t))), t G [0,r]. On account of (gjZTj) , 
([222D, and (I3TT2|1 we thus find 

±( d n Q+ \(t v 0) - f B(t) ( driQ+ )(t v 0) 

d-J (0 '^ 0) = U 

where A(t) = H'J p (p(t)) and B(t) := B(t, y*) is defined as in (ET23|) with y = y+, 
for t G [0, r]. On the other hand, 

<»°> i(aM»))=(.i _%)w«»(a(Mo)), 

Hence, using = d x (V x H(x,V(p)) = H" x (x,V(p) + H" p (x, V<£>) <p" on if and 
(16. 1 01) fc (16. 11 1) , we can verify directly that 



X(t)\ / dpX(t,p(t)) 
Pit)) ' V-^"(7(t))d p X(t,p(t))+d p P(t, 



ie [0,- 



solves 

Jt\p)"'~ \ » -B(/i T ; \p)"' \p)" n ' \i 
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which implies (16. 9p . □ 
Lemma 6.4. The following identity holds, 
H , -uj ^ _ d A 0*,2/*) d_1 det (exp^exp" 1 ^*))) 



v x * id v Q + (r, yir ,0)J \V(x*)\\V(y*)\(l -V 2 (x*))^(l -V 2 (y+))^' 

Proof. We drop the subscripts * of the distinguished points x*, in this 
proof. First, we introduce some notation. Let bi,...,bd denote some G(y)- 
orthonormal basis of R d such that b d — (1 — ^ /2 (2/)) _1 ^ 2 ^y/|^| and ci, . . . , q 
be some G(x)-orthonormal basis of R d such that q = (1 — V 2 (x))~ 1 ' 2 v x /\v x \. 
Let 6*, . . . , b d and c*, . . . , c d denote the corresponding dual bases and let B 
and C denote the matrices whose z-th row is b* and c*, respectively. Then we 
have Bv y = (1 — V 2 (y)) 1 ^ 2 \v y \ e d and C v x = (1 — V 2 (x)) 1 ^ 2 \v x \ e d , where e d is 
the <i-th canonical basis vector of R d . Therefore, we find, using (I6.9p . that is, 
id v Q + (r, y, 0) = d p X(r, y, w(0)), 

det (J c) det (v x <Q^y,0)) det (J 
= (1 - ^(x)) 1 / 2 H (1 - ^(y)) 1 / 2 K| det ( ° _ . „ 



Since det 5 = (1 - F 2 (y)) d / 2 , detC = (1 - V 2 (x)) d ' 2 , 
and, analogously, = (1 — y 2 (s)) 1 / 2 /|l / (a;)|, we obtain 



(6.12) 



(| , -uj" ^ _ det((c*rf p X(r,y, W (0))(6*) T ) 1 ^ K(i _ 1 ) 



w x id v Q+{r,y,0)J \V(x)\\V(y)\(l-V 2 (x))^(l-V 2 (y))^' 

In order to compare the position space projection, X, of the Hamiltonian flow 
associated with H and the exponential map at y we observe that 

X(r,y,p) = exp, [d A {X(r,y,p),y) (1 - V^y))' 1 ' 2 p/\p\ 

for p e R d in some neighborhood of £7(0), since £(p) := (1 — V 2 ^))^ 1 ^ 2 p/\p\ is 
normalized with respect to G(y) and the initial momentum p of a Hamiltonian 
trajectory is collinear with its initial velocity in our case. We set r(p) := 
d A (X{r,y,p),y). Then it follows that 

d p X(r, y,p) = exp' y (r(p) £(p)) [£(p) <g> r'(p)] + r(p) exp' y (r(p) £(p)) £'(p) . 
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By GauB' lemma we know that c* exp' y (r(w(0)) £(tu(0))) £(ro(0)) = 0, for i = 
1, . . . , d — 1, thus 

c* d p X(r, y, 07(0)) (&*) T = r(w(0)) c* expj, (r(«7(0)) £(w(0))) C'(w(0)) (b*) T , 

for j = 1, . . . , d— 1. If -P^( ) denotes the orthogonal projection in R d onto the 
Euclidean orthogonal complement of ro(0), we have 

C\m(0)) (6*) T = (1 - F 2 (y))- 1/2 |^y| P^ (0 ) (1 - V 2 (y)) 6, = 6, , 

since |w(0)| = (1 — V 2 {y)) 1 ^ 2 . Using also the identities r(zu(0)) = dx(x,y) and 
r(ro(0)) £(oj(0)) = exp~ 1 (x) we arrive at 

det ((c*d p X(r,y,w(0))(b*) T ) 

= d A {x,y) d ~ l det ((c* exp^exp" 1 ^)) ■ 
Finally, we use c^exp^exp" 1 ^)) bd — 1 to conclude that 

det ((c* d p X(r,y,o7(0)) (6J) T )i<ij^rf-i) = d A (x, y)^ 1 det (exp^exp" 1 ^))) . 
Inserting this identity into (1 6 . 1 2 j) we arrive at the assertion. □ 

Appendix A. Connection to the BMT equation for Thomas 

precession 

In this appendix we consider only the case d = 3 and choose the standard 
representation of the Dirac matrices, 

We know that the matrix-valued term M(x+,y+) = U(r) (— V(y*)) A + (izu(0)) 
appearing in f l 1 . 1 3 j) maps RanA + (itu(0)) onto RanA + (itjj(r)). If we choose 
appropriate bases of these subspaces then the coefficient matrix of M(x*,y+) 
corresponding to these basis vectors is a solution of some spin transport equa- 
tion which is closely related to the Bargmann-Michel-Telegdi (BMT) equation 
for the Thomas precession of a classical three-dimensional spin; see, e.g., [UE]. 
In our special case, where no magnetic field is present, the Hamiltonian deter- 
mining the particle trajectory Qj is H given by (HI, and y/1 - w 2 = -V(y), 
this spin transport equation reads 

d o- ■ (E(x) x p) 

(A.l) - B (t) = im(7(t)Mt))s(t), M{x,p):=- 



dt v ' ww ' w/ v / ' ■ -2V{x)[l-V{x))Y 

where E = — W is the electrical field and s(t) is a complex (2 x 2)-matrix. 
Notice that the usual momentum is replaced by an imaginary momentum in 
H since we are dealing with some sort of tunneling regime. The BMT equa- 
tion, or rather its analogue for the Hamiltonian H, is obtained from (lA.lj) by 
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choosing some u G C 2 and computing the differential equation satisfied by the 
expectation value s(t) := (s(t)u\(rs(t)u) C 2 of the vector of Pauli matrices. 
In our case the BMT equation turns out to be 

d s(t) x (Ejrm) x zu(t)) 



In order to derive flA.lj) and connect it to (11.131) we start with Equation (I4.28P 



for B+, which, on the domain @ + , reads 

(A.3) - 2V(x) A + (W<p(x)) (id t + cx ■ V x ) B° + (t, x, 0) = , (t, x, 0) G ®+ . 

For every x G R 3 , the range of A + (iVy3(x)) is spanned by the two mutually 
orthonormal eigenvectors of D(x,iVip(x)) in the (4 x 2)-matrix 

(A 4) W(x) = 1 (I 1 ~ V ^)\ 1^ 

y/-2V(x)[l-V(x)] V a ■ zV ^( x ) 

so that A+(iVp(a;)) = W(x) W(x) T . Next, we write B%(t, x, 0) = W(x) C(t, x), 
for some (2 x 4)-matrix C(t,x). This is possible since fulfills (To). Then 
(1A.3|) implies IU T (id t + ck • V x ) W C = and we observe as in pQ that the 
operator —iW T (idt + ot ■ V x ) W equals 

d t + F-V x + ^divF-iVJl(x, V<p) , F(x) := V p H(x, V<p(x)) = V(x) -1 Vp(2;) . 
Substituting 7 for x and using 7 = ^(7), we deduce that 

|c(t, 7 (t)) = ~divF( 7 (t))C(t l7 W) + ia«(7(t),w(t))C'(t.7(t))- 

Notice that, unlike (I4.35p . the previous equation involves the complete diver- 
gence of F. Using the ansatz C(t,j(t)) = g(t) s(t) W(y±) T , where g is scalar, 
and 5 solves (lA.lj) with 5(0) = I2, we thus find by means of Liouville's formula 
and the initial condition (7(0,?/*) = W(yi,) T that g(t) = det[d y Q + (t, y+, 0)] -1 / 2 , 
t G [0, t}. We arrive at the following formula for B+ alternative to (14.341) . 

£?°(r,x*,0) = det[d y Q + (t,y^)Y l,2 W(x,)s(T)W(y,y. 

Using the previous formula instead of (14.341) in the proof of Proposition 16.11 we 
obtain the following asymptotics for the Green kernel, 

_ (1 - U 2 (x,))i(l - V 2 (y^ (VjxJVjy*))* (1 + Q(h)) e - d ^^' h 
det [exp^exp- 1 ^*))] 172 2n d A (x*,y+) /h 

■ W(x+) s(r) W{y*) T , s solves fl£3), 5(0) = 1 2 . 
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Here W is given by (IA.4j) where Vyj(:c*) = tu(t), V(p(y±) = zu(0). We remark 
without proof that the scalar term in the second line is the asymptotic ex- 
pansion of the Green kernel of the Weyl quantization of a/1 + £ 2 + V in three 
dimensions. This can be inferred by means of a procedure similar to the one 
carried through in the present paper. Finally, we remark that s(r) can be rep- 
resented by means of the polar coordinates of a suitable solution of (IA.2I) and 
additional dynamical and geometric phases; see [U §4]. 

Appendix B. Proof of Lemma 13.31 

In order to give a self-contained construction of the phase functions ip± we 
present the proof of Lemma I3T31 in this appendix. The proofs below are variants 
of those in [7j where the symbol is assumed to be homogeneous of degree one. 
We recall the definition of & in (13.81) and start with the following lemma which 
corresponds to [7J Lemma 1.7]. 

Lemma B.l. For every compact subset K C Q, we find some Ck G (0, oo) 
such that, for all (s,p) = (s,x,£) G C x K, 

(B.l) (£a ± e)( S ,p) > -^a ± (Mp)-C K \%p\ 3 . 

Proof. Let a be a + or a_. Since J^ a is a real differential operator we have 
[J#a, 3?] = 0, whence 

X(s+ (x|9fcf )) = a- (V € o|0 + (V s a|9?f) - (x\m x a) 

= a — ( V^a | i3£ ) — ( i^sx | V ' x a ) — K( x \ V ' x a ) 

on Q. Taking the imaginary part and using that [J^, 3] = we obtain 

— Jf^ & = 3 (a — ( V p a | i^sp )) . 

Taylor expanding both a and V p a at Kp using d^^^VpafJRp) = 0, a G Ng d , 
we infer that 

(B.2) a-(V p a\z%p) = a(Kp) + \ ( 3p | a'^p) 3p ) + 0(\%p\ 3 ) . 
Furthermore, a Taylor expansion of 3a(3ftp ± tSp) with t > yields 
i Sa(Kp) + l - ( 3p | 9aJ p (3fip) 9p> 

(B.3) = — (3a(Kp + t3p) + 3a(Kp-t3p)) +tO(|Sp| 3 ) . 

The 0-symbols in flB.21) and flB.31) are uniform when p varies in a compact 
subset of Q. Choosing t — 2 and using that 3a ^ we thus arrive at the 
assertion. □ 
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Proof of Lemma \3.3[ We drop all sub- and superscripts ± in this proof. Taylor 
expanding the right side of 4 ^sK t = ^sJt^i^t) at 9Rft f and using Duhamel's 
formula we obtain, exactly as in pp. 351], the estimate 

(B.4) \%k u \ < 0(l)(\QK t \ + I \^a'(UK r )\dA . 

J u 

It holds for all p in a compact, complex neighborhood of (yo,r]o) and u,t G 
[0, r+£i], for some £i > 0. Since, for p in a compact set, the curves [0, r+£i] 3 
t i — y ^R,K, t {p) stay in a compact set, we may apply the standard estimate for 
positive functions to flB.4p . which together with Holder's inequality gives 

(B.5) < 0(1) (IQfKtl + (4) 1 / 2 ) , 4(p) := / -%a(MK v (p))dv, 



for ^ u ^ r ^ t ^ r + £x- Next, we integrate the estimate (IB. ip for 
^(3(Q,/tt) = J^((5)(q, Kf) from w to t, use f lB.5j) to bound |9te r |, r e [w, £], 
and arrive at 

(B.6) 6( fi , «*) ^ 6(^, /c) + ^ - 0(1) (t - «) (l^l 3 + (/*) 3/2 )- 
In the case u = we further observe by means of (13.71) and flB.5j) that 
(B.7) 6( ft , «o) ^ -O(l) |^p| 3 ^ -O(l) (|^| 3 + (J*) 3 / 2 ) . 

If r > in the plus-case we use the assumption ^sa + (y , rj ) = and Lemma 
to deduce that 3a + (re+(y , rj )) — 0, t e [0, r], thus I*($/o,7?o) = 0, O^w^t^ 
r. Back in the general case we conclude that, for ^ w ^ t ^ r + £j, we can 
ensure that the integrals I l u {p) are as small as we please by assuming that p 
is contained in a sufficiently small neighborhood of {yo,f]o), and that e± > is 
sufficiently small. Then ( 1B.6j) and ( 1B.7j) give (13. 9p and the estimate 

(B.8) 6(q, «*) ^ 6( fe , k u ) 0{l) {t - u) \%K t \\ 

for 0^u^t^T + Ei. Squaring (1B.5|) with r = u, dividing by some suitable 
constant, and adding the result to (IB. 81) we obtain 



©(ft, «t) + -\%n t \ 2 > 6(<r„, ««) + ^k u \ 2 + ^- - -J /* - 0(1) (t - u)|3^| 3 , 

where we can assume that the constant (coming from (IB. 5^ ) satisfies C ^ 
4. By possibly restricting the neighborhood of (yo,fjo) further to ensure that 
maXf e [o )T + £l ] \5sKt \ is sufficiently small, we infer from the previous estimate that 

(B.9) 6(ft, K t ) + |S^| 2 ^ ©fo, k u ) + i |3«; u | 2 + ~ 4 , 

for 0^M^t^r + £i. On account of (13. 9ft . where the integral is non- negative, 
we may again restrict the neighborhood of (yo, rjo) so that |3ffc u (p) | is sufficiently 
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small to ensure that (1 — ^)&{<; u , k u ) + ^\^sn u \ 2 ^ 0, for all u G [0,r + £i\. 
Subtracting the latter term from the right hand side of flB.9j) and using C ^ 4 
we arrive at 

(©(q, «*) + \%n t \ 2 ) > 1 (©(^, « u ) + |3/g 2 ) + i /* , 

which is f)3.10p . Subtracting -^(&(q r , n r ) + | | | 2 ) ^ from the right hand 
side of (13.101) . where the integral is non-negative, we finally obtain (13.111) . □ 
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